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ABSTRACT 

A  coma-corrected  zoned  mirror  has  previously  been  shown  by  para- 
geometrical  optics  to  posess  no  spherical  aberration  and  no  third  order 
coma.   However,  geometrically  defined  aberrations  provide  no  information 
about  the  diffraction  pattern  structure  which  may  be  of  importance 
in  microwave  applications.   In  this  investigation  a  study  using  a 
more  rigorous  boundary  value  approach  has  been  made  of  a  zoned  mirror 
with  emphasis  on  its  effectiveness  in  coma  correction.   Because  a 
cylindrical  mirror  has  no  astigmatism  and  other  higher  types  of  aber- 
ration (according  to  Nijboer' s  classification)  and  because  the  mathemati- 
cal analysis  is  simpler,  a  two  dimensional  zoned  mirror  was  chosen  for 
analysis. 

The  image  patterns  (the  field  intensity  distributions  on  the  focal 
plane  for  various  incident  angles  as  commonly  used  in  optics)  as  well 
as  the  main  features  of  the  radiation  pattern  (minor  lobe  level,  gain 
and  beam  width  for  various  scan  angles  as  commonly  used  in  antenna 
engineering)  are  obtained.   These  features  are  compared  with  those  of  a 
smooth  parabolic  cylinder  with  same  focal  length  and  aperture.   It  is 
found  that  for  a  system  with  small  F  number  and  nearly  uniform  illumination, 
the  zoned  mirror  is  much  superior  to  the  smooth  parabola  in  coma  correction. 

The  gain,  minor  lobe  level  and  beam  width  of  the  radiation  pattern  for  a 

o 

zoned  mirror  are  relatively  constant  for  a  scan  angle  up  to  25  ,  while 

those  for  the  corresponding  smooth  parabola  change  considerably.   The 
ratio  of  the  two  secondary  maxima  of  the  intensity  distribution  in  the 
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focal  plane  (a  function  of  the  coma  aberration  only)  is  substantially 
equal  to  unity  in  the  case  of  a  zoned  mirror  whereas  it  changes  drasti- 
cally for  a  smooth  parabola. 

Calculations  of  the  diffracted  fields  are  based  upon  three  different 
orders  of  approximation  for  the  induced  current  on  the  mirror,  namely, 
(i)   Geometrical  optics  current  distribution   (ii)  a  distribution  taking 
into  account  the  effect  of  the  edges   (iii)   a  distribution  taking  into 
account  both  the  edge  effect  and  the  first  order  coupling  between  zones. 

It  is  found  that  the  field  intensity  on  the  focal  rlane  (and  by 
reciprocity,  the  far  field)  obtained  by  using  (ii)  and  (iii)  are  not  signifi- 
cantly different.   Results  obtained  by  using  approximation  (i)  follow 
those  obtained  by  (ii)  and  (iii)  quite  closely.   Thus  the  simple  and 
straight-forward  physical  and  optics  solution  seems  to  offer  a  result 
of  good  accuracy  even  with  a  large  number  of  discontinuities  in  the 
diffracting  surface  as  in  the  present  case.   However  for  higher  accuracy 
the  edge  effect  should  be  considered. 
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1.   INTRODUCTION 

For  wide  angle  scanning,  the  aperture  of  a  microwave  reflector  is 
always  limited  due  to  the  progressive  deterioration  of  the  image  as  the 
aperture  of  the  system  is  increased.   Although  the  gain  and  directivity 
increase,  the  scanning  performance  deteriorates  because  of  spherical 
aberration,  coma^  astigmatism  etc.  which  are  inherent  with  wide  aperture 
systems. 

The  correction  of  optical  errors  for  wide  aperture  systems  is 
generally  performed  by  an  assembly  of  lenses  which  provides  the  designer, 
for  each  element,  with  two  radii  of  curvature  and  the  index  of  refraction 
as  the  three  available  parameters  for  correction  of  optical  errors.   In 
the  case  of  a  reflector^  however^  only  the  radius  of  curvature  is  available 
as  a  parameter  for  correction  of  optical  errors.   The  correction  of  aberra- 
tions by  placing  the  surfaces  of  one  curvature  on  the  surface  of  another, 
such  as  a  zoned  mirror,,  has  been  suggested  by  various  workers.   This 
method  is  of  practical  interest  and  importance  due  to  its  simpler  con- 
struction and  lower  cost  as  compared  with  a  lens  system  particularly  in 

VHF  and  UHF  ranges. 

28 
Ramsay  and  Jackson   have  made  an  experimental  investigation  of  a 

zoned  mirror  particularly  in  comparison  with  a  spherical  mirror  in  scan- 

27 
ning  performance.   Lately  Toraldo  di  Francia   and  his  coworkers  have 

shown  by  his  generalized  geometrical  optics  approach  that  coma  aberration 

is  greatly  reduced  in  a  zoned  mirror.   The  purpose  of  this  investigation 

is  to  analyze  the  performance  of  the  zoned  mirror  as  a  wide  angle  scanning 

antenna  by  employing  a  more  rigorous  boundary  value  problem  approach. 


2.  COMA- CORRECTED  ZONED  MIRROR 

2,1   Description  of  the  Zoned  Mirror 

Let  us  consider  a  set  of  confocal  paraboloids  of  focal  Lengths 
differing  by  X/2  having  the  common  focus  at  c  and  axis  VC.  (Figure  1). 
Then  the  family  of  the  paraboloids  will  be  defined  by 

2(f  -  n  h 

r=- —  ,     n  =  0,  1  2,    ...  (1.1) 

1  +  cos  41 

Where  f  is  the  focal  length  of  the  paraboloid  with  largest  focal 
length. 

This  family  of  paraboloids  will  be  equiphase  in  so  far  as  all  rays 
parallel  to  the  axis  VC}    after  reflection  from  all  the  parabolic  surfaces, 
will  meet  at  the  focus  c  in  phase. 

These  confocal  paraboloids  will  intersect  the  spherical  surface  2, 
having  c  as  the  center  and  f  as  the  radius^  as  circles  represented  by  P, ,  P„ 
etc.   If  a  set  of  planes  perpendicular  to  the  axis  VC  is  drawn  through  the 
points  N  ,  N   ....  representing  the  intersection  of  the  paraboloids  with 
the  axis  VC,  they  will  cut  the  spherical  surface  2  in  circles  represented 
by  Q  ,  Q   „...  etc,   A  family  of  coaxial  cylinders  with  VC  as  the  axis  and 
cross  sections  given  by  the  circles  Q  ,  Q   etc.  are  drawn.   Then  the  parabo- 
loidal  surfaces  cut  out  by  these  circles  form  a  zoned  mirror. 

Thus  A, B1 ?    A  B_  etc.  are  the  reflecting  zon^s  made  of  perfect  conductor 

while  EL  A„ ,  B„A„  etc.  are  the  sections  of  coaxial  cylinders  marking  the 
1  2'   2  3 

boundaries  of  A  B  ,  A  B  etc.  which  are  assumed  to  play  no  role  in  image 

formation. 

It  can  be  proved  that  all  the  tangents  to  these  zones  at  P,  .  P  etc, 
^  &  12 


FIGURE    1 


4 
pass  through  V.   Thus  the  parabolic  zoned  surfaces  can  be  approximated  by 
portions  of  a  set  of  cones  having  the  same  vertex  V. 

2C2   Some  Earlier  Works  on  a  Zoned  Mirror 

28 
Ramsay  and  Jackson    investigated  experimentally  the  performance  of  a 

coma  corrected  zoned  mirror  similar  to  the  one  described  above  with  the 

modification  that  the  tips  and  not  the  centers  of  the  zones  were  on  the 

spherical  surface  (coma  sphere) .   This  was  done  presumably  for  the  purpose 

of  simplifying  the  construction  and  it  was  assumed  by  the  authors  that 

this  difference  did  not  change  the  performance  to  a  great  extent.   Ramsay 

and  Jackson  compared  the  measured  gains,  side  lobe  levels  and  beam  widths 

of  the  radiation  pattern  of  a  zoned  mirror  with  those  of  a  spherical  mirror 

o 
with  a  feed  offset  of  15  .   They  found  that  the  zoned  mirror  is  slightly 

better  than  the  spherical  mirror  but  this  advantage  is  doubtful  when  the 

difficulty  of  construction  is  taken  into  consideration.   However^  it  is 

felt  that  they  worked  with  a  case  where  the  advantages  of  a  zoned  mirror 

are  greatly  suppressed  due  to  the  following: 

(1)  The  F  number,  namely  f/D  ratio  used  by  them  is  unity,  where 

f  is  the  focal  length  and  D  is  the  diameter  of  the  aperture  of  the  assembly. 
Thus  the  dishes  used  for  comparison  by  them  are  rather  shallow, 

(2)  A  15  db  tapered  illumination  is  used  by  them.   This,  in  effect, 
weights  heavily  a  small  central  portion  of  the  dish  which  is  essentially 
the  same  for  both  the  zoned  mirror  and  the  spherical  mirror, 

(3)  The  zones  are  not  arranged  to  make  their  mean  surface  coincide 
with  the  so  called  "coma  circle.''   In  fact  except  for  the  first  zone,  none 
of  the  zones  has  any  point  lying  on  the  coma  circle. 


(4)  A  three  dimensional  zoned  mirror  as  used  by  these  authors,  if 
coma  corrected,  will  have  astigmatism  as  the  main  aberration.  In  their 
experiment,  they  seem  to  have  made  no  effort  in  differentiating  these 

aberrations . 

27 
Ronchi  and  Toraldo  di  Francia   by  employing  the  theory  of  "Parageo- 

metrical  Optics,"  *  have  obtained  expressions  of  the  fifth  order  coma  and 

astigmatism  for  this  type  of  mirror.  (The  third  order  coma  is  absent  as 

a  result  of  the  satisfaction  of  the  sine  condition.)   Toraldo  di  Francia, 

very  recently  reported  that  their  experimental  results  are  in  good 

agreement  with  the  theoretical  prediction  of  the  best  focal  surface  for 

scanning  purpose. 


*   Briefly  speaking  the  parageometrical  optics  is  essentially  the  theory  of 
geometrical  optics  generalized  for  the  grating  systems  which  ascertains 
the  relation  of  the  diffracted  ray  vector  in  terms  of  the  incident  ray 
vector  and  the  grating  system  geometry. 


3.   OUTLINE  OF  THE  PRESENT  INVESTIGATION 

3.1   Limitations  of  the  Aberration  Study  by  Optical  Methods 

Aberrations  of  any  system  are  generally  studied  by  either  (1)  geometri- 
cal optics,  namely  the  ray  tracing  technique  or  (2)  obtaining  the  diffraction 

pattern  at  the  image  space  with  an  assumed  aberration  function  at  the 

25 
aperture.     The  drawback  of  (2)  is  obvious  in  that  the  over-all  perform- 
ance of  a  given  system  is  still  unknown.   About  (1),  it  is  well  known  that 
geometrical  optics  loses  its  effectiveness  at  microwave  frequencies  since 
there  exis  ts  no  definite  boundary  between  the  shadow  and  the  illuminated 
regions,  i.e.  the  image  is  no  longer  a  well  defined  geometric  one.   Thus, 
aberration  defined  in  terms  of  geometrical  optics  loses  its  clearness  of 
description  at  lower  frequencies.   This  is  particularly  true  for  the  region 
near  nulls  and  minor  lobes  where  even  the  physical  optics  method  begins  to 
show  errors. 

It  is  difficult  to  translate  the  geometrical  optics  evaluation  of 
aberration  into  a  diffraction  pattern.   However,  since  the  structure  of 
minor  lobes  and  coma  effect  are  closely  related  (the  ratio  of  the  two 
minor  lobe  levels  on  either  side  of  the  major  lobe  is  a  function  of  the 
conventionally  defined  coma  aberration),  an  investigation  of  the  minor 
lobes  seems  to  be  in  order. 

For  a  system  like  the  zoned  mirror,  having  a  large  number  of  edges 
and  discontinuities  in  structure,  the  physical  optics  method  seems  to  need 
further  justification.   It  is  reasonable  to  assume  that  the  diffraction 
due  to  the  edges  and  coupling  between  zones  could  be  important  particu- 
larly in  the  regions  of  nulls  and  minor  lobes  upon  which  the  coma  aberration 


has  a  dominant  effect. 

3.2   Plan  of  the  Proposed  Investigation 

In  recent  years  due  to  the  rapid  growth  of  microwave  technology,  great 
interest  has  been  developed  in  the  generalized  diffraction  theory  and  many 
workers,  e„g.  Fock  *  Braunbeck  ,  Keller      etc.  have  propounded  various 
theories  for  approximate  solutions  of  general  diffraction  problems.   It, 
therefore,  seems  that  a  complete  system  of  relatively  simple  structure 
can  be  analyzed  with  reasonable  accuracy  by  the  "boundary  value  problem" 
approach. 

The  surface  of  revolution  described  in  2.1  presents  severe  mathematical 
difficulties,  however.   Therefore  a  two  dimensional  version  of  the  zoned 
mirror  with  each  zone  approximated  by  a  strip   is  chosen  for  this  investiga- 
tion.  In  fact  this  system  is  of  more  interest  since,  being  free  from 
astigmatism,  it  enables  us  to  evaluate  its  ability  in  coma  correction 
exclusively, 

It  is  therefore  proposed  to 

(1)  Analyse  the  cylindrical  system  using  the  "Boundary  value  problem" 
approach  so  that  more  detailed  information  regarding  the  performance  of 
the  system  can  be  obtained  than  is  possible  by  other  qualitative  methods 
mentioned  before,, 

(2)  Provide  a  comparative  study  of  possible  alternative  approximations 
in  the  analysis. 

(3)  Provide  information  on  antenna  performance,  such  as  gain,  side 
lobe  level,  beam  width,  for  various  scan  angles. 


*   The  approximation  of  each  zone  by  a  strip  is  valid  for  all  zones  except 
the  central  one,  for  which  a  correction  has  been  used  in  later  sections 


3.3   Formulation  of  the  Problem 

The  problem  can  be  stated  as  follows:   Given  (a)  a  two  dimensional 

zoned  mirror  consisting  of  a  number  of  perfectly  conducting  plane  strips 

(1  2}     ...  n) ,  infinitely  long  in  z-direction  arranged  as  described  before 

(Figure  2)  and  (b)   a  plane  electromagnetic  field  E     with  the  electric 

inc 

field  polarized  parallel  to  the  edges  of  the  strips^  i.e.  having  only  the 
z  component,  incident  on  the  zoned  mirror;  it  is  required  to  find  the 
intensity  distribution  of  the  electric  field  on  the  image  plane. 


Y 


n 


FIGURE  2 


When  a  plane  wave  is  incident  on  the  zoned  mirror  each  element  of  the 
zoned  mirror  will  have  a  current  induced  on  the  surface,  the  distribution 
of  which  will  depend  on  the  angle  of  incidence  and  the  width  of  the  strips 
It  will  also  depend  on  the  current  distribution  on  the  surface  of  other 
strips  because  of  the  coupling  between  the  strip  under  consideration  and 
other  strips,   Once  the  current  distribution  on  each  strip  is  determined 
the  electric  field  at  any  point  in  space  can  be  evaluated. 

This  problem  can  not  be  solved  by  the  well  known  method  of  solution 


9 
of  the  wave  equation  by  separation  of  variables  since  the  assembly  can  not 
be  defined  by  a  complete  coordinate  surface. 

However,  it  is  possible  to  have  an  integral  equation  formulation  of 
the  problem.   Boundary  value  problems  of  a  conducting  or  dielectric  body 
in  a  given  electromagnetic  field  may  be  formulated  by  expressing  the 
scattered  field  in  terms  of  an  integral  of  induced  sources  over  the  surface 
or  throughout  the  volume  of  the  scatterer.   The  boundary  conditions  at  every 
point  on  the  scatterer  then  lead  to  an  integral  equation  in  terms  of  the 
induced  source  function.   The  integral  equation  formulation,  unlike  that  of 
differential  equations,  contains  the  boundary  condition  specified  in  the 
formulation,  but  generally  is  not  traceable.   It  does  not  require  the  system 
to  be  separable  and  often  leads  to  an  approximate  solution  more  directly,, 

Let  the  currents  induced  on  the  surface  of  the  zones  1,  2,  3  ...  n 
etc.,  be  denoted  by  I  ,  I   ...  I   etc.,  respectively. 

Then  the  vector  potential  has  only  z  component.   A  at  any  point  P 

S5 

due  to  these  currents  I.  ,  I  ....  I   etc.   is  given  by 

12       n 

r 

A  a         I  (s)g(s,s')ds  +   !      ln(s)g(s,s/  )ds  4-  ...  +        I  (s)g(s,s)ds 
z    J        1  j       2  J      n 

strip  1  strip  2  strip  n 

(3.1) 

where  s(x,y)  is  a  point  on  the  strip  under  consideration,  s  (x,y)  is  the 

point  of  observation,  and  g,  the  free  space  Green  's  function  is  defined  by 


g  *  |-  H  (2)  Ik  /(x-x')2  +  (y~y')2  k  with  k  =  B.         (3.2) 
4i   o        '  A. 

From  Maxwell    s   equation   we   have, 

H=VXA  (3.3) 
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and 


Hence 


which  in  our  case  gives 


E=^r    VXH 
ioo€ 


E   -  -±-  V  X  VX  A 

ico€ 


A  1  A 

E      z    -  — -VX  V  X  A      z 
z  icot  z 


(3.4) 


(3.5) 


(3.6) 


or 


E      z   =  -    (VV-    A      z    -    z  72A    )    =  ^—  A      z 
z  y  z  z  icot      z 


(3.7) 


Since 


and 


2  2 

V      A     +    k      A     =   0 
z  z 


■3-  A      =   ° 
OZ       z 


At  the  surface  of  any  strip,  the  total  tangential  electric  field  must 
be  zero.   Hence  the  diffracted  tangential  field  at  the  surface  must  be 
given  by 


E   «  =  E 
Diff    z 


(3.8) 


inc 


Therefore 


^—   .  A  s - E 
i<^€    z      z 


(3.9) 


inc 


at  the  surface  of  any  strip.   Thus 


k 

ico€ 


I    (s)g(s,s')ds 
strip   1 


I2(s)g(s;,  s    )ds   + 


I    (s)g(s    s    )ds 
n 


strip   2 


=    -    E 


(3.10) 


Jinc 
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The  above  represents  n  equations  for  s   lying  on  any  of  the  n  strips. 

Representing  the  operator 

2 

g(s,  p)ds 


ioo€ 


by  G   we  can  express  Equation  (3C10)  symbolically  as 

qp 


GI  =  -  E 


(3.11) 


where 


G  = 


Gll   G21 


Jnl 


G     G    .  .  .  „  G 
12     22        n2 


In     2n 


nn 


I  = 


Equation  (3.11)  represents  a  set  of  n  simultaneous  integral  equations  which 
when  solved  gives  the  values  of  the  current  distributions  I,  ,  I_,  ....  I 
on  all  the  n  strips. 

Let  I.  i,e.; current  on  any  strip  i  be  represented  as 


I.  =  1. .  +   2  I. 
1     ll    j=l   1J 


(3.12) 
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where  I    is  the  current  induced  on  the  surface  of  the  strip  i  due  to  an 
11 

incident  plane  wave  in  the  absence  of  other  strips,  and  I. .  is  the  current 
induced  on  the  surface  of  the  strip  i  due  to  the  field  of  I .  of  the  strip  j 

The  solution  of  Equation  (3.11)  is  very  difficult  to  obtain^  but  a 
method  of  successive  approximation  can  be  used  to  obtain  the  solution  to 
any  desired  accuracy,  at  least  in  principle. 

Choosing  the  i    equation  from  the  set  denoted  by  the  Equation  (3.11) 
we  have  for  s   on  i 


icoi 


n 
I  (s)  (s,s  )ds  +   S 

1  j=l    - 

strip  l  strip  j 

j^i 


I  (s)  g(s^s/)ds 


=  -  E 


inc 


(3,13) 


substituting  Equation  (3.12)  in  Equation  (3.13)  we  get 


n 

s 
j=l 


since 


strip  i 


I .  .  (s)  g(s,  s  )ds  + 


I . (s)  g(s, s  )ds 
J 


strip  j 


for  s'  on  i 


=  0 


(3.14) 


ioo£ 

strip  i 


I.  .  (s)  g(s,s')ds  =  -  E 


(3,15) 


which  defines  the  integral  equation  for  the  isolated  strip  i.   Using 
Equation  (3.12)  again  in  Equation  (3.14)  we  obtain 


N 
j=l 


r  / 

I   (s)g(s,  s  )ds 
strip  i 


+    f  I  .  (s)g(s,s')ds  +  S         I 

J  .        .    JJ  k=l   J  .   .  " 


:  (s)g(s  ,)ds 


strip  j 


Mj 


strip  j 


(3,16) 


when  s   is  on  i 
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where  I  .  is  the  current  induced  on  strip  i  due  to  strip  j,  (j  ^  i)  and 

I.,  is  the  current  induced  on  strip  j  due  to  a  strip  k^  (k  ^  j).   Since 

s   is  on  i,  g  in  the  first  integral  of  (3.16)  is  much  greater  than  that 

in  second  and  third  integral  and  further  I  .  >>  I  ,  .   Thus  an  iterative 

J  J     jk 

procedure  can  be  established  by  solving 


n 

S 

j=l 


I^)(s)g(s,s/)ds  +    J      I   (s)g(s,s')ds 


strip  i  strip  j 


=  0 


/  (3.17) 

for  s   on  l 


Equation  (3,17)  represents  an  independent  linear  integral  equation  and 
there  are  in  total  n  similar  equations  obtainable  by  putting  i  si,  2, 
3,  ....  n. 

Thus  the  set  of  n  simultaneous  integral  equations  represented  by 
Equation  (3.11)  have  been  reduced  to  n  independent  equations  for  each 
order  of  the  solution^  a  typical  one  of  which  is  given  by  Equation  (3.17). 

Since  Equation  (3.17)  is  linear  and  I  .    represents  the  current 

induced  in  strip  i  due  to  I . .  we  get 

JJ 

I.J(s)g(s>s/)ds  +    I       Ijj(s)g(s,s/)ds  =  0 

strip  i  strip  j 

(3.18) 
for  all  j  ^  i 

Equation  (3018)  can  be  solved  if  I    is  known.   In  that  case  it  represents 

JJ 

the  integral  equation  of  the  current  distribution  on  a  strip  i  due  to  a 

known  source  distribution  I .  .  on  strip  j  „   Further  by  definition  I  .  .  is 

JJ  JJ 

the  current  induced  on  strip  j  due  to  a  plane  wave  incident  on  it.   Hence 
Equation  (3.15)  represents  the  integral  equation  for  I... 


14 

Thus  the  problem  has  been  reduced  to  the  solution  of  Equation  (3.15) 

and  using  the  value  of  I  ..  there  from^  to  solve  Equation  '3.18).   Equation 

(3018)  can  be  solved  if  the  two  dimensional  Greens  function  i.e.  the  current 

I    on  strip  i  due  to  a  filament  source  of  unit  strength  is  known, 
ij 

Thus  Equation  (3.15)  is  a  special  case  of  Equation  (3.18).   The 
rigorous  solution  of  both  of  these  equations  leads  to  a  series  solution 
involving  Mathieu  functions.   There  are  also  various  approximate  solutions 
of  sufficient  accuracy  for  the  two  above  equations  and  these  have  been  dis- 
cussed in  the  next  chapter. 

(1) 

Once  I.  .  and  I.  .   have  been  evaluated  by  solving  the  above  equations, 

the  current  distribution  on  all  the  strips  can  be  determined  and  putting 

the  values  of  I.(s)  in 

l 

2    n     r 
-r—     Z  I       I.(s)g(s,s')ds  =  E(s)  (3.19) 

iW€i=l   J 

strip  i 

the  field  at  any  point  s  can  be  calculated. 
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4,   BRIEF  SURVEY  OF  PREVIOUS  WORK  ON  THE  DIFFRACTION  BY  A  STRIP  OR  A  SLIT 

This  has  been  the  object  of  investigation  by  several  authors. 
There  are  at  least  two  methods  which  lead  to  the  exact  solution,  viz 

(1)  Method  employing  separable  coordinates  to  solve  the  wave 
equation. 

(2)  Integral  equation  method. 

4.1   Method  Employing  Separable  Coordinate  System 

20 
Morse  and  Rubenstein  obtained  the  solution  of  the  problem  of 

diffraction  by  a  strip  due  to  an  incident  plane  wave  using  the  separable 

coordinate  systems.   The  result  is  obtained  in  terms  of  infinite 

series  of  Mathieu  Functions.   The  summary  of  their  work  is  given 

below. 

Let  a  perfectly  conducting  elliptic  cylinder  be  considered  with 

its  axis  parallel  to  the  electric  vector  E.    of  an  incident  plane 

mc 

wave  whose  direction  of  propagation  makes  an  angle  a  with  the  major 

iwt 
axis  of  the  ellipse.   In  all  discussions  the  time  factor  e    will  not 

be  written  but  understood. 

Let  the  distance  between  the  two  focii  of  the  ellipse  be  2h„ 

Using  an  f-lliptic  cylindrical  system  of  coordinates,  we  have. 


x  ~  h  cosh  fe  cos  t|   y  s  h  sinh  \    sin  n  and 


z  ~   z 


and  the  diffracting  ellipse  is  defined  bv  &  .   The  wave  equation 

o 
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FIGURE  3 


expressed  in  an  elliptic  cylindrical  coordinate  system  is  given  by 


Q   P       Q   P  p   O 

— ~o   +   — o   +  2  k  h  (cosh  2^  -  cos  2T| )  E  =  0 


(4.1) 


Assuming  a  solution  of  the  above  in  the  form  <j>(^)i|j  (H )}     (4„1)  is  reduced 
to 


E± 


+  (s  cosh  b-b)  <j>  -»  0 


(4.2) 


and 


— 5"  +  (b-s  cos  rj)  i|j  _  o 

an2 


(4,3) 
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2  2 
where  b  is  the  separation  constant  and  s  =  k  h  .   The  solutions  of 

(4.2)  and  (4.3)  are  given  by 


oo 


Se  (s,ix)  =  2  Dev  '        cosh  (2k+p)x 
r        k=0   2k+p 


and 


OO 


(r) 


So  (s,ix)  =  i   2   Do„,    sinh  (2k+p)n 

k=o   2k+p 


(even  solution) 


(4.4) 


(odd  solution) 


(r) 
Se  (s,x)  =   2  De„,    cos  (2k+p)x 

r  }  2  k+  p 

k=0      p 


(even  solution) 


(4.5) 


So  (s  x)  =  2  Do*    sin  (2k+p)n 

k=o   2k+p 


(odd  solution) 


where  p  =  0  if  r  is  even   and    p 
and 


1  if  r  is  odd 


OO 


Se  (s,0)  =  1     i.e. 

r  ' 


2  De     =1 

k=0   2k+p 


using  the  boundary  condition  i.e.  E^    on  the  surface  of  the  diffracting 

tan 

ellipse  £>      must  vanish,  Morse  and  Rubenstien  obtained  the  diffracted 

o  ' 

field  given  by 


18 


ED 


iff 


inc 


Je„_(£    ) 


27T     S 


••  (2tl 


n=0   I  "2n  He""t^    ) 
2n      o 


.    He;    ;(e>-Seo    01)Seo    (a) 
2n  2n  2n 


.-L-^f^Ho^d)   Sog        «>    So,  fa, 

N2n+2  Ho0(2)(l    )  2n  2n+2  2n+2 

2n        o 


+   J 


No      t    t,    (2)     /f    \         2n+l  2n+l  2n+l 

2n+l   He„      ,  (b    ) 


2n+l      o 


(4.6) 


+  ^h  7^44  Ho^(e>  s°2-i01)  s°2-i(a) 

n+1   Ho  (§    ) 

2n+l      o 


where  En .  f  f  is  the  diffracted  field  at  any  point  P^,1!)  and  E.    is  the 

incident  field.   When  £   =0,  the  elliptic  cylinder  degenerates  into  a 

strip  of  width  2h  and  Equation  (4.6)  gives  the  scattered  field  due  to 
a  strip. 

By  using  the  relations    curl  E  =  -  B    and    I  =  n  X  H 


where  n  is  the  outward  normal  unit  vector  to  the  strip,  the  current 


density  at  any  point  on  the  surface  of  the  strip  is  found  to  be 
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4Er 


I    = 


inc 


'77 


-       Se2n(T1>    Se2n(a> 
2Li 


jcyh  sintl  |   yT    1n=0  (2) 

L  2n  2n 


+    J 


5e2n+l(T1)    Se2n+l(a) 

N        ,    He;    '       (0) 
2n+l         2n+l 


+      2 
n=0 


rSo  CH)    So„      9(a)  So  Cn)    So0         (a) 

2n+2  2n+2  2n+l 2  n+1 

(2)  ;  :        +  J  „/       „  (2) 


N  n      „  Ho;-'  (o) 
2n+2    2n+2v 


NQ   .  Hox"'  .  (o) 
2n+l     2n+l 


(4.7) 


If  the  incident  field  instead  of  being  plane  is  that  due  to  a 
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filament  parallel  to  the  z  axis  with  strength  I   located  at  a  point 

S(£  j1!  ),  it  can  be  shown  that  the  Hertz  potential  at  any  point  P^,1!) 
o  o 

is  given  by 


7T  = 


—  S~ 
-iu>£  n=0 


hI2)(§  )  ze  (£)  Se  01)  Se  01  ) 
n    o     n      n      no 

TT) 

N  He  '  (o) 
n    n 


H*2fe 


(4.8) 


&J 


n   o 

+   /  (2) 
N  Ho  (o) 
n   n 


zo  (£)  So  01)  So  (11  ) 

n      n      no 


for  the  region  defined  by  0  <  §  <  £   .   For  the  region  defined  by  5  ?  5 

—   —  o  o 

interchange  g  and  4   in  Equation  (4.8). 

o 

Further 


ze  (u)  =  Ne  (u, )  Je  (u)  -  Je  (u  )  Ne  (u) 
n        n  1    n        n   1    n 


(4.9) 


zo  (u)  =  No  (u, )  Jo  (u)  -  Jo  (u  )  No  (u) 
n        n  1    n        n   1    n 


(4.10) 


n 


je  (u)  =  rn_ 


-u 


2(De  ) 
o 


°°      l~Z  khe 

2  rib  (-X)     D62r  V"^ > 


(4.11) 


if  n  is  even 


Je  (u)  = 
n 


'77 


r-- 


n-1 


-u 


u 


-u 


2(De  )    r=0 


2De      rj  (Sll)j     (i^I)  +  j     <»*•) 

2r+l   r    2     r+1   2        r+1    2 


,khe  .     .  .        ,, 

j  ( )     if  n  is  odd 

r   2 


(4.12) 
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r-5  ,,    u         .,    -u 


Ne    (u)    =  /77  «f  ,      .        2^        v    ,khe    .  _     .khe 

/ J        S  (_1)  De2r   ri~2~)Jr(~~2~~)  n    1S    GVen 

2 (De    )         r=o 


n-i 
2  „  ,  ,     -u  ,  ,    u  , ,     -u  , ,    u 


=    Fjf~  2  (-1)  De        .  r,    .khe      v  .khe    .  .khe      .        .khe    , 

— o  2r+l[J    ( — - )Y      _  (-^ )+J      ,( — = — )Y„(-^ — ) 

2      r=o  r        2  r+1      2  r+1        2  r      2 


2(De1) 


if  n    is    odd  (4.13) 


n 
r— =•  ,  ,    ~u  ,  ,    u  ,  ,    ~u  u 


l/  2(Doz)2      r=l  r      r_1        2  r+1      2  r+1      2  r"1      2 

if  n   is   even 

"""T"  r.    .khe"U.  .        ,kheu,.        .khe"",.    ,kheu 


if  n   is   odd  (4.14) 


v2      r=o 
2(DQl) 


No    =  rn_  2  (-D 


n 
•f   ,      .r~2~  r  khe"uv  ,kheuv  ,khe~Uv  ,kheU 


2r      r-1        2  r+1        2  r+1        2  r      2         J 


2(D)2      ™ 

o  .  _ 

if  n   is   even 

n-1 
r  — - u  u  -u  u 

iZ ZZ  ,   ,  ^         2   _  r  _    .khe      .  ,khe   .     _         ,khe      .        .khe    .  . 

It s(-D         D,2r+1[jr(_T_)Yr+i(_r)-j    i(___)1fr(-T.)] 

2(Dl)2      r=° 

if   n   is   odd  (4.15) 


He      (u)   =   Je    (u)    -    i   Ne    (u)  (4.16) 

n  n  n 


Ho2)(u)    =   Jo    (u)    -    i    No    (u)  (4.17) 

n  n  n 


N     =  7T[2    (De    )      +    (De    )2   +    .  ..]         if  n   is   even 
n  L       v      o  2 


(4.18) 


OS  Q 

=  V  2     (Deo      _)  if  n   is   odd, 

2  r+1 
r=o 
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From  Equation  (4.8)  and  (4.9)  the  current  induced  on  the  surface  of  the 
strip  can  be  obtained  as 

_     h!2)(£    )Se    (n)Se    01    )         H&?^    )Se    01) So    fl ) 


I  =  -I       2 
o 


n        o        n  no 


lVVaen 


(2) 

n=o  N      He    (o) 

n        n 


N'   Ho^o) 
n        n 


(4.19) 


4.2   Integral  Equation  Method 
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Schwarzchild    solved  the  diffraction  problem  due  to  a  slit  in  a 

plane  infinite  screen  which  is  complementary  to  a  strip.   His  results^ 

2 
obtained  in  a  different  way  by  Baker  and  Copson  .  are  indicated  below. 

Let  the  screen  and  slit  occupy  the  plane  x  =  0  in  rectangular 

cartesian  coordinates.   The  edges  of  the  slit  of  width  2h  coincide  with 

lines  x  =  0}    y  =  0  and  x  =  0,  y  =  -2h.   Let  the  incident  electric  vector 

E.    be  polarized  parallel  to  the  z  axis  and  the  normal  to  the  wave  front, 

inc  } 

i,e.  direction  of  propagation  make  an  angle  a  with  the  positive  x-axis. 
Then  H   is  the  only  nonvanishing  component  of  the  magnetic  field  tangential 

y 

to  the  screen  and  a  current  of  density  I  (y,a)  flows  on  the  screen  in 

the  direction  of  the  z  axis.   The  total  E  field  at  any  point  (x  ,  y  )  is 

o   o 

then  given  by 


kz 
„  ,     .     ik(x^  cos  a+y   sin  a)     o 

E  (x  y  )  =  e     o       'q       — 

zoo  4 


-2h 


1  ,  nxH  <2)(kr)dv 
z(y,a)  o 

(4.20) 


where 


/xq2  +  (y-yo)2 


and 


z  =   I  u. 
O   ■'  C 


When  the  boundary  condition,  i.e.  E  (x   y  )  is  zero  at  all  points  on  the 

'      '  zoo 
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screen,  is  applied  the  following  integral  equation  is  obtained. 

-iky  sin  a 
4e    ° 


+  I    hz(y'a)  Ho(2)(k!y-yoi)dy  = 


kz 


(4.21) 


for  y  >  0  and  y  <  -2h. 
o  o 

Using  Fox's  Formula   and  changing  the  variables,  Equation  (4.21)  can 

be  transformed  into  a  pair  of  simultaneous  integral  equations  of  Wiener- 

Hopf  type  having  solutions  given  by 


4 
and 


L*  I  (y  a)  *  d,Ay  *>  -  I  f  !H  I  (-y-2h  a)  e  -ik(y+z+2h>  /^2h 


dz 


(4.22) 


o  T  ,   „,   .      -i2kh  sin  a         1 
—  Iz(-y-2h,a)  =  e  i|i(y,-a)-  - 


z  g-ik(y+z+2h) 

'  I  (v  a)  — /z+2h.  .„ 

>a;     y+z+2h   /(-— — )dz 


4    z 


where 


,7T 


,  „   .    cos  a-   iky  sin  a.   e 
4*(y,a)  *  — - —  e  -  cos 


iky  sin  a    r—- —   .    7T   a. 
a  e  y  F[  v^2ky  sin  (-+-)] 


if 

-:L4    -iky 
e     e   '      .V        a. 

+ sin  ( -) 

4   2 


(4.23) 


y/2ff   /ky 


where  F(u)  is  the  Fresnel  integral. 

Equations  (4.22)  are  exact  and  can  be  solved  by  successive  approxi- 
mations which  yield  a  solution  of  the  form 


I  (y,a)  =  S  ln  (y,a) 


(4.24) 


n=0 
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where 


,ftO 


1   f  -ik(y+z+2h)    — 

I   (y,a)  =  -  =      I  -,  (z,a)  S /z+£n  dz    for  n  >  i 

nz  ^   j    n-lz   »       y+z+2h      V  y 

^°  (4.25) 

and 

I  (y,a)  =  i|>(y,a). 
o 
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The  series  represented  by  Equation  (4.24)  has  been  proved  by  Schwarzchild 

to  be  uniformly  and  absolutely  convergent  for  all  values  of  h  greater  than 

zero  if  k  is  real. 

Schwarzchild  evaluated  the  first  term  of  the  series  of  Equation 

,         29 
(4.24)  which  physically  corresponds  to  considering  Sommerfeld  s  formula 

for  two  half  planes  considered  separately  for  the  slit  and  neglecting 

the  diffraction  by  each  half  plane  of  the  field  scattered  by  the  other. 

Solutions  of  the  same  problem  in  series  form  in  increasing  power  of  the 

29 

ratio  of  slit  width  to  wavelength  have  been  obtained  by  Sommerfeld, 

and  Bouwkamp.    These  solutions  are  particularly  effective  when  the  slit 
is  relatively  narrow. 

The  results  obtainable  from  the  solution  in  Mathieu  functions  have 
their  use  restricted  to  the  availability  of  appropriate  Mathieu  function 
tables.   Further  these  functions  converge  very  slowly  far  large  arguments. 
Thus  when  the  strip  width  is  small  conpared  to  wavelength  the  Mathieu 
function  solutions  are  very  useful.   The  same  is  true  for  Sommerfeld  s 
method  in  which  the  rate  of  convergence  of  the  series  depends  on  the 
narrowness  of  the  strip  or  slit. 


24 

30 
The  integral  equation  method  employed  by  Schwarzchild  on  the  other 

hand  is  good  for  very  wide  strips  in  which  case  the  convergence  is 
extremely  fast.   In  fact  with  strips  or  slits,  large  compared  to  wave- 
length, even  the  first  term  of  the  series  as  obtained  by  Schwarzchild 
corresponding  to  the  first  order  approximation  gives  very  satisfactory 
results. 

For  slits  and  strips  not  large  compared  to  the  wavelength,  Karp 

13  19 

and  Russek   and  Millar   have  obtained  approximate  solutions.   Their 

solutions  are  very  suitable  when  the  width  of  the  slit  is  neither  narrow 

enough  for  Sommerfeld's  series  solution  nor  wide  enough  to  use  the  first 

term  of  Schwarzchild ' s  series  solution. 

13 
4,2.1   Approximate  Solution  by  Karp  and  Russek 

The  solution  obtained  by  Karp  and  Russek  has,  as  the  first  term,  the 
term  calculated  by  Schwarzchild  from  the  exact  series  solution  and  a 
correction  term  which  accounts  for  the  interaction  between  the  two  edges 
of  the  half  planes  forming  the  slit.   The  effect  of  the  mutual  interaction 
between  the  two  edges  are  accounted  for  by  lumping  these  effects  in  the 
form  of  two  fictitious  filament  sources  at  the  two  edges  of  the  slit  and 
then  determining  the  strength  of  these  two  filaments. 

The  effectiveness  of  the  solution  arises  from  the  fact  that  the 
strength  of  the  fictitious  line  sources  at  the  edges  of  the  slit  which 
are  used  for  representing  the  effect  of  diffraction  due  to  the  half 
planes  can  be  obtained  in  terms  of  simple  trigonometric  functions  and 
Fresnel  integrals, 
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4„2„2  Asymptotic  Solution  by  Millar 

H.  F„  Millar  obtained  an  asymptotic  representation  for  each  term 

of  the  series  solution  of  the  integral,  equations  given  by  Equation  (4„22) 

in  terms  of  the  inverse  powers  of  kh,   He  has  shown  that  each  term  of  the 

series  of  Equation  (4„24)  gives  rise  to  an  asymptotic  expansion  whose 

-1 
leading  term  is  of  the  order  (2kh)  2  compared  to  that  of  the  preceding 

term.   The  asymptotic  expansion  of  the  series  solution  to  any  finite  order 

in  inverse  powers  of  (2kh)  can  therefore  be  obtained  by  considering  a 

finite  number  of  terms  of  the  series, 

Millar  has  calculated  a  few  higher  order  terms  of  the  asymptotic 

expansion.   The  calculation  becomes  very  involved  as  the  order  of  term 

is  larger0   He  has  calculated  the  transmission  coefficients  taking 

different  number  of  terms  of  the  asymptotic  series  and  comparing  them 

31 
to  those  obtained  by  Skavlem„     A  study  of  this  table  reveals  that  for 

some  moderate  values  of  2kb  (e  g0  10  i.e   of  width  of  about  1„6X)  the 

best  result  is  obtained  when  only  the  first  term  of  the  expansion  is 

considered „   It,  therefore,  seems  reasonable  to  expect  a  satisfactory 

result  even  if  only  the  first  term  of  the  series  solution  is  used  for 

widths  larger  or  equal  to  one  wavelength, 

21 
4„2„3  Approximation  by  Moullin 

Moullin  and  Phillips   have  computed  the  current  distribution  on 

strips  of  different  widths  for  incident  plane  waves  at  different  angles 

from  the  Mathieu  function  solution  given  by  Equation  (4,7)  and  compared 

numerically  to  that  near  the  edge  of  a  half  plane,,  for  the  same  incident 

plane  wave.   They  found  that  the  difference  between  the  two  becomes  very 
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small  for  strips  over  X/2  wideu   Therefore  they  concluded  that  the 
current  distribution  near  the  edge  of  a  strip  is  sensibly  independent 
of  the  presence  of  the  other  edge  i.e.  the  coupling  between  the  edges 
is  very  small  if  the  strip  is  at  least  one  half  wavelength  wide,,   Hence 
they  approximated  the  current  distribution  over  an  isolated  strip  by  that 

which  would  appear  if  the  strip  were  a  part  of  an  infinite  plane  with 

2  2 

two  filament  sources  located  at  the  two  edges  of  the  strip,   Moullin 

also  reached  the  same  conclusion  for  a  cylindrical  wave  diffracted  by  a 
strip. 
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5„   APPROXIMATION  USED  FOR  THE  SOLUTION  OF  THE 
DIFFRACTION  PROBLEM  OF  A  ZONED  MIRROR 

It  has  already  been  shown  that  the  density  of  current  induced  on  the 
surface  of  a  ribbon  due  to  a  plane  wave  incident  at  an  angle  a.  with  the 
E  vector  parallel  to  the  edge  of  the  ribbon,  in  elliptic  cylindrical 
coordinates  defined  by  b  and  "H  is  given  exactly  by 


1  K^lnfo  K,<  (0)     '   N    Hef>   (0) 

2n  2n  2n+l   2n+l 


(5.1) 


?   Se2n+2(11)So2n+2(a)     .  S°2n+1  (T1)S°2n+l  ^ 

n=o    N     Eov    J       (0)         N    Hov  '   (0) 
2n+2   2n+2  v  '  iN2n+l   2n+l  K    ' 


2  2 

S  =  k  h  }    2h  being  the  width  of  the  ribbon  and  r\      determining  the 

location  of  a  point  on  the  ribbon,, 

In  order  to  determine  1 01) ,  the  Mathieu  functions  Se  ,  He   and  the 

'  n'    n 

constants  N  have  to  be  computed  for  each  of  the  ribbon  widths  2h.   This 

2  2 
can  be  easily  accomplished  when  s  i0e„  k  h   is  small,  since  the  Fourier 

expansion  of  the  Mathieu  function  series  converge  rapidly  under  that 
conditionQ   However  when  s  is  large,  (corresponding  to  a  strip  a  few 
wavelengths  wide)  as  is  the  case  for  some  strips  in  the  present  problem, 
the  usefulness  of  the  Mathieu  function  series  is  lost  because  the  converg- 
ence of  the  series  is  extremely  slow.   The  other  difficulties  involved  in 
the  computation  of  this  function  with  large s  have  been  discussed  in  the 
appendix.   It  has  been  shown  there  that  the  computation  of  current  density 
across  the  surface  of  a  wide  ribbon  using  Mathieu  function  solution 


28 

Equation  (5.1)  is  not  only  involved  and  impractical,  the  accuracy  obtainable 
is  not  the  same  for  all  values  of  *\    and }       in  fact,  it  is  quite  low  when 
T|  approaches  7T/2. 

On  the  other  hand  a  physical  argument  indicates  that  as  the  strip 
becomes  very  wide,  the  solution  for  each  half  of  the  strip  should  approach 
that  of  a  half  plane  which  is  readily  obtainable  in  terms  of  Fresnel 
Integrals. 

Millar,  Karp  and  Russek,  and  Moullin  and  Phillips  in  one  way  or 
another  obtained  solutions  along  this  approach.   They  extablished  their 
theory  by  demonstrating  numerically  the  small  difference  between  the 
approximate  and  exact  solutions  for  a  certain  set  of  parameters  for  which 
the  latter  is  available,  without  going  into  a  general  treatment  of  the 
error  involved  in  the  approximation.   In  fact  the  first  order  term  of  the 
solution  obtained  by  all  of  them  is  the  same  as  that  obtained  by  Schwarz- 
child  which  is  sufficiently  accurate  for  strips  of  the  order  of  a  wave- 
length wide. 

5,1   Approximation  for  the  Current  Distribution  across  an  Isolated  Strip 
Due  to  an  Incident  Plane  Wave 

As  discussed  in  the  last  section,  for  a  strip  of  width  no  less  than 

a  wavelength,  the  first  order  Schwarzchild  Js  solution,  namely  the  Sommer- 

feld  s  half  plane  solution  furnishes  a  satisfactory  answer  to  the  problem. 

The  current  distribution  across  a  strip  can  therefore,  for  the  present 

purpose,  be  approximated  in  the   following  way. 

A B 

B 


FIGURE  4 
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Let  the  strip  AB  have  a  current  distribution  across  its  surface  given  by 
I  (Figure  4) .   We  then  assume 


I  = 


( 


Iv,  +  I>,  -  I       on  the  strip 


I 


0  outside  the  strip 

where  In  is  the  current  distribution  across  a  half  plane  extending  from 

A  to  oc  . 

1^   is  that  current  distribution  across  a  half  plane  extended  from 
B 

B  to  -to  and  I   is  the  current  distribution  across  an  infinite  plane. 

OO 

All  these  current  distributions  Iv,  ,  IL   and  I   are  known  exactly 

UA'   "B      oo 

in  terms  of  Fresnel  Integrals  and  simple  functions.   Thus  I  can  be  easily 
computed  if  this  approximation  is  used.   With  this  approximation  it  is 
noted  that  current  across  the  half  plane  with  edge  at  B  and  that  across 
the  infinite  plane  are  nearly  equal  at  points  far  off  from  the  edge  B, 
Hence  if  AB  is  large,  the  current  distribution  near  the  edge  A  of  the 
strip  will  approximately  be  that  due  to  the  half  plane  with  its  edge 
at  A0   The  same  is  also  true  near  the  edge  B  of  the  strip  AB.   While  at 
the  middle  of  the  strip,  since  lu  as  Iu  =s  I    the  current  is  approximately 
I  ,   This  agrees  closely  with  the  conclusions  of  Moullin  and  Phillips. 

The  error  involved  in  the  above  approximation  with  regard  to  matching 
the  boundary  condition  at  the  surface  of  the  strip  has  been  calculated 
and  found  to  be  very  small.   Further  the  computation  from  the  asymptotic 
expansions  of  the  exact  solution  as  discussed  in  the  appendix  are  both 
involved  and  inaccurate.   Thus  the  above  approximation  probably  does  not 
yield  results  inferior  to  that  one  could  obtain  from  the  exact  solution 


by  asymptotic  expansion. 


Evaluation  of  I 
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p(r,0) 


FIGURE  5 

-    ,,       .   .,   .   ,  ,         ikr  cos  (<b-<j>  )     .     ,/N 

Consider  an  incident  plane  wave  given  by  E  =  e         ^  Y   =  u  *r,  cb;d>  ) 

z  o 

with  its  E  vector  parallel  to  the  edge  of  the  half  plane  defined  by  the 
lower  half  of  y-z  plane.  r}ty   are  the  coordinates  of  the  point  of  obser- 
vation  P}    and  <}>   is  the  direction  of  the  incident  wave.,  all  measured  from 
the  origin  A  as  shown  in  Figure  5,   The  total  E  field  is  given  by  Soramer- 
f eld ' s  half  plane  solution  viz 


E  (r. §;§'   )  =  u(r;<j);(j)/) 


-  u(r,<j);3  7r-<])) 


(5.3) 


where 


*4  r  T 

e  I    -  -i  \  9 


v/¥'   o 


and 


J 


T  =  /2kF  cos  ^ 


2 
(5.4) 
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Therefore 
j  IT 


|/? 


ikr  cos  (<(>-<j)/  ) 


/2kr  cos  *^x- 


i-iX2(.^  _  eikr  cos  (<|)-37r+<f>) 


-00 


V2kr  cos  ^ — _i£ 


e    d\ 


(5.5) 


If  the  angles  are  measured  from  the  half  plane,  then 

IT  / TT  / 

a  =  —  +  <j>    and    a  =  -s   +  <J> 


and 


4 


,/2kr  cos 


a -a 


e   ,  ikr  cos  (a -a  ) 

E  = [  e 

VI  J  CO 

■*'  — oo 


-iX2  -     ikr  cos  (a+a  ) 
e    dX  -  e 


./Ski- 


cos 


a+a 


—  X  K,2t  ,  <* 
e    d\ 


-oo 


(5.6) 


Using  the  relation  H 


iW|j. 


curl  E,  the  magnetic  field  is  obtained„   On 


the  surface  of  the  strip,  the  surface  current  density  induced  due  to  the 

incident  electromagnetic  wave  is  given  by  the  discontinuity  of  the 

magnetic  field  i.e. 

i   =   x   X  [H         ,         -   H         ,  ] 

tan        n        tan        0_ 
a  ~0  I  a=27T 


H 


1  z 


y  7 


tan        iwu,   3x 


9e 


=  2xy      ~-  U-f 

Ua\x      ox 


JhA 


a=0 


9E 
z 

3x~| 


a-2fl" 
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from    (505) 


3E 
z 


e       1 


/? 


kr    ,    ,     a   -a  a+aN    -ikr 

(sin  -y-  -   sin  —  )e 


/     ikr  cos  (a -a)         ,       /s 
+   ikr  [  e  'sinCa-a) 


-iX2   ,  -ikr   cos(a+a)    .  / 

e  dX  -   e  sin(a+a) 


'-oo 


'e-lX2dX 


where 


?1   -     /2k7 


cos 


a-a 


and 


a+a 


/ 


T      =      /2kr   cos 
2  2 


It    therefore   follows    that 


Ti 
i- 

2e        I  /     ikr   cos  a 

Iv,  .   -  \  k  sin  a   e 

nA 


. Q. 

j/2kr   cos  — 


l/ffwfi 


i       , a 

-   \/2kr   cos  - 


-iX2  .         .    /2k  a7    -ikr     ~  ^ 

6  /       dX  "    Vr"  Sin  2   e  I   XXy 


(5.7) 


I,    ,    which  flows   parallel   to  the  edge  of   the  half  plane  has   its   phase 
referred    to  A0      Now  consider 

// 2kr   cos  —  ^/2kr   cos  — 

VlXadX  .  2    J  e-iX2dX 

-  \J2kr   cos  —  -"o 

Substituting  X2  =  |  t2  we  get 


r/2krcos|  2  y^  cos  |       ^ 

e"iX  dX  =     pa 


J 


e'V  dt  =  fa  [cA-isA] 

V  2 
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where 


C.    = 


7T      2 

cos  —  t   at 
2 


and 


SA  = 


IT      2 
sin  —  t  dt 
2 


and 


o      Ar  a- 

u  =   2    /—cos  - 


Using  the  above  we  get 

2 
h^  =Tp-" 


I-      =  — [  sm  a    ^  (C  +S    )cos(kr   cos  a)-(C   -S    )sm(kr   cos  a)  \   +  sin  — 

I  J        /TTkr 


(cos   kr   -   sin  kr) 


+  — —  [sin  a    i  (C  +S    )sin(kr   cos  a)+(C   -S    )cos(kr   cos  a)f   - 


1  a 

sin  — 

fkr 


(cos   kr  +   sin  kr) ]  (5.8) 

1^      is   obtained   by   substituting    (7T-a)    for  a      and    r   for   r      in  Equation 
(5.8),    since    (r?    Tl-a.)    is    the   observation   point   referred    to  B.      Thus 


Iv^  =  —  [sina   J  (C  +S    )cos(kr   cos  a)+ (C   -S    )sin(kr   cosa)i  + 
LIB       T|  IBB  B      B  J 


1      cos  a 


'7Tkr 


(cos   kr   -   sin  kr) ] 


+  —  [sin  a       (Cd-S  .)    cos  (kr   cos  a)-(C  +S   )sin(kr   cos  a)      -  cos  — 


|/7Tkr 
(sin  kr  +   cos   kr) ] 


(5.9) 


where 


CB  =     I        cos  |   t2   dt,       SB  =  sin  |   t2   dt    and   V  =   2    /p      sin  | 
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The  phase  of  Iy.   is  referred  to  the  point  B.   I   can  be  calculated 

nB  oc 

in  a  similar  way.   It  turns  out  that 

2k   .    /  ikr  cos  a  ,_    . 

I   =  —  sin  a  e  (5.10) 

where  r  is  the  distance  of  the  point  of  observation  from  the  edge  A  and 
the  phase  of  I   is  referred  to  A.   Thus  it  has  been  possible  to  evaluate 

OO 

in  terms  of  Fresnel  Integrals  and  simple  trigonometric  functions, 

I(=I^  +  IL  -  I  ),  the  current  distribution  at  any  point  on  the  surface 

of  a  strip  due  to  a  plane  wave,  incident  at  an  angle  —  -  a.   This  expression 

does  not  take  into  account  the  proximity  of  other  strips.   Thus  the 

current  I  is  what  has  been  called  L  . ,  the  self  induced  current  of  the 

strip  i,  previously. 

5.2   Evaluation  of  the  Mutually  Induced  Current 

The  effect  of  the  incident  wave  on  an  assembly  of  strips  is  to 
induce  a  self  current  in  each  of  the  strips.   The  proximity  of  the  strips 
will  influence  the  current  distribution  on  the  other  strips,  i.e.  there 
will  be  a  current  induced  in  each  strip  due  to  the  scattered  field  of 
other  strips,  which  may  be  called  mutually  induced  current.    Thus 
the  resultant  current  distribution  in  each  strip  is  the  sum  of  the  self 

tl       v 

and  mutual  currents. 

The  mutually  induced  current  can  be  found  by  successive   iteration 

a  ii 

as  discussed  previously.   In  order  to  evaluate  the  mutual   current  on 

any  strip  due  to  another  strip  in  its  neighborhood  consider  a  strip  AB 

and  a  filament  source  of  strength  I  located  at  S.   The  current  induced 

on  AB  due  to  the  source  at  S  can  be  evaluated  exactly  in  terms  of  infinite 
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summation  of  products  of  Mathieu  functions  as  given  by  Equation  (4.19) 
of  the  previous  chapter.   However  the  limitation  of  this  solution 
discussed  there  applies  to  this  case  as  well.   We,  therefore,  seek  an 
approximate  solution  suitable  for  numerical  computation.   Following 
the  case  of   self  current,   the  current  distribution  across  AB  can 
also  be  approximated  by  the  sum  of  the  current  distributions  on  the 
two  half  planes  with  edges  at  A  and  B  less  by  the  current  distribution  on 
a  doubly  infinite  plane  each  under  the  influence  of  the  filament  source 
at  S,  as  before. 


y 


^ 


FIGURE  6 


The  part  due  to  infinite  plane  can  be  obtained  simply  by  the  image 
method.   The  part  due  to  the  half  plane  can  be  obtained  in  two  alternative 
forms:  the  series  solution  involving  Bessel  functions  and  the  integral 
expression.   Now  consider  the  second  form  first.   We  know  that  the  field 
at  any  point  P  due  to  a  source  of  strength  I  at  S  in  presence  of  a  half 
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plane  with  edge  at  A  (Figure  6)  is  given  by 


18 


iQJu.1 


r  o 


z  47T 


-ikR   cosh  £   c 
e  db    - 


A 


-ikR,  cosh  £  ,c 
el  db 


-OO 


J 


-OO 


(5.11) 


where 


0-0 


Sinh  £„  -  2    /rr      cos  — - 
0  y o  2 

R 


and 


e+e] 

Sinh  £      =  2  ,/rT  cos  —— 
1  /      o  2 

R, 


R  =7r     +    r        -   2rr      cos  (0-0   )      and      R     -Jr     +   r     - 
o  o  o  IV 


2rr      cos (0+0   ) 
o  o  o 


From  the   relation   -zH  =   curl   E     where  z   =   i^u.,    the   tangential    component 

of   the  magnetic   field   due   to  E      is 

z 


tan 


6e 
i  z 

o3jI  730 


Now 


TO  d£ 

9  -ikR   cosh  |   *  -ikR   cosh  £n        0 


0 


30 


•d£ 


j         30 

■OO 


-ikR  cosh  § 


d| 


-7 o   (r+r   )    sin 


0-9        -ik(r+r   ) 


R 
ikrr 
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r  • 


sin    (0  -   0   ) 
R  o 


cosh 


t      ikR   cosh  £  ,t 
be  db 
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..    H 


Be 

i  z 


tan       co\xr       9© 


I_ 
4tt 


[r~     (r+r   )  0-0        -ik(r+r   )      ikrn 

-°-.   °     sin  °  e  °  -  _°.  sin    (0-0  ) 

V  r  „2  2  R  o' 


R 


-00 


,     c      -ikR      cosh  £   t 
cosh  be  d£ 


/  r      (r+r   )  0+0        -ik(r+r    )         ikr 


-ikR^    cosh  | 
cosh  |   e  d£ 


Hence 


H 


tan 


0=0 


I_ 

27T 


rr^    (r+r   )  0        -ik(r+r   )      ikr 

o  o'  o  o'  o      .      A 

—  sin  —  e  -  — - —  sin  0 

2  2  R  o 


2<Tc 

'       R 
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-oo 


0 
COS    — 
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-ikR  |/z  +1    . 
e  '  dz 


and 


H 


tan 


0=27T 


I_ 

27T 


r      (r+r   )  0        -ik(r+r   )         ikr 

o  o  _ 

-  sin  0 

R  o 


o  o        .        o 

sin  — —  e 


-2/^  0 

— r cos  2 


Vr  R2 


-00 


-ikR  /z2+l    . 
e  dz 

by  change  of  variable  sinh  §  =  z. 

The  density  of  the  current  induced  on  the  surface  of  the  half  plane  with 

edge  A  due  to  E   is  given  by 
z 


/     A 

IhA  =  x  X 


tan    tan 
o      27T- 
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I12  =  j  l'dl  (5.16) 

CD 


where  di  is  an  elementary  length  of  CD  and  I  is  the  Green  s  function 

for  the  induced  current  at  a  typical  point  on  strip  (1)  due  to  the  current 

at  another  typical  point  on  strip  (2)  as  given  above. 

The  alternate  form  of  the  solution  may  be  obtained  for  the  induced 
current  1^  and  1^  by  method  of  separable  variables  in  cylindrical 

A  D 
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coordinates   as  discussed  in  Chapter  3.   They  are 

XhA  s  h    (A+1  •>  (5'17) 


ih'B  ■  h'^+i  2)  (5°18) 

where 

A  =  S  (-1)  (2n+l)Y   !(27rr)J   .  (27TR)  cos  — ±_  a               (5.19) 

n=0            n+2      n+2  2    o 


and 


C=  2  (-l)n(2n+l)J   !(27Tr)J   ,  (27TR)  cos  ^J^  a  (5.20) 

n^O  n+5      n+2 


where  a   =  7T  -  8    when  R  >   r  interchange  r  and  R  in  A  and  C. 
o        o 

Similarly  for  In<   k'  and  C  are  obtained  by  replacing  r,    R  by  r  and  R 
respectively  where  r'and  r'  refer  to  B  as  the  origin  i.e.  for  the  half 
plane  extending  from  B  to  -oo„   With  this  formula,  one  of  the  double 

integrations  required  by  the  other  method  is  avoided  by  taking  summation 

"  ii 

instead.   Either  of  the  two  will  give  the  value  of  the  mutual  current 

correct  to  the  first  order. 
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5.3   Investigation  of  the  Error  Due  to  Assumed  Current  Distribution 

From  the  discussion  in  the  previous  chapter  it  is  clear  that  with 
the  assumed  current  distribution  across  a  strip  the  scattered  field  at 
any  point  will  not  be  exactly  the  same  as  would  be  obtained  from  the 
exact  current  distribution.   It  is  therefore  necessary  to  estimate  the 
change  caused  by  the  assumed  distribution  over  the  exact  one  so  far  as 
the  scattered  field  is  concerned.   In  order  to  estimate  this  we  need  to 
know  the  field  produced  by  the  exact  current  distribution  which,  in  turn, 
implies  that  we  must  know  the  exact  current  distribution.   However,  in 
this  section  an  expression  has  been  obtained  by  which  it  is  possible  to 
estimate  the  effect  of  the  assumed  current  distribution  on  the  tangential 
E  field  over  the  surface  of  the  strip.   The  total  tangential  E  field  on 
the  surface  of  the  strip  with  the  exact  current  distribution  must  be 
zero.   The  total  tangential  E  field  on  the  surface  of  the  strip  with 
assumed  current  distribution  will  therefore  provide  an  estimate  of  the 
closeness  of  the  assumed  current  distribution  to  the  exact  one0 

Let  the  exact  current  distribution  at  any  point  x  on  the  strip  AB 
be  given  by  I  (x)  „   Then  the  field  at  any  point  P(r,  0)  due  to  the  current 
distribution  I  (x)  on  the  strip  AB  is  given  by 

2    rh 

^   \         I  (x)H  ^[kfr-xlldx^  E  (r,0)  (5,21) 

46t)£   J     so1'    lJ      z  > 

-h 

When  the  point  of  observation  P  is  on  the  surface  of  the  strip  AB,  the 

the  total  tangential  E  field  on  the  surface  being  zero,  we  have  the 

integral  equation  for  I   given  by 

s 
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~-  I   I  (x')  H*2*(k|x-x'l)  dx'  =  -  E,   (x)  for  -h  <  x  <  h.    (5.22) 

4««  Is  o  inc  —   — 

-fa 


The  assumed  current  distribution  I(xO  will,  in  general,  not  satisfy 
the  above  integral  equation.   Suppose  I(xO  satisfies  a  slightly  different 
integral  equation 


h 

|   I(xO  H(2)(klx-x'l)  dx  =  -  [E,   (x)  +  AE(x)] 
l         o  inc 


-b 


for  -  h  <  x  <  h        (5.23) 

where  A  E(x)  is  the  amount  by  which  the  right  hand  side  of  (5.23)  differs 
from  that  of  Equation  (5.22).  Combining  Equations  (5.22)  and  (5.23)  we  get 

fa 

~T2  f   (I  (xO  -  I(x')>  H(2)(klx-x"i)  dx  =  AE. 

4Ufc   /      s  o 

-h 

Thus  A  E  indicates  the  error  committed  in  the  boundary  condition  by  taking 
the  assumed  current  distribution  I(xO  for  the  problem  of  the  strip. 

Although  I  (xO,  the  exact  current  distribution  is  an  unknown  function, 

s  ' 

the  left  hand  quantity  can  still  be  found  due  to  the  fact  that  I(xO  con- 
sists of  half  plane  and  whole  plane  solutions  which  satisfy  certain  equa- 
tions. 

As  stated  before  IhA(xO,  IhR^x  "^ >    and  ■"■  (x')  are  solutions  of  half 
plane  and  whole  plane  problems,  and  they  satisfy  the  following  equations. 

h 

-■Z'~£  K    (x">  H*2*[klx-x'l]  dx<"  «  -  E„       for  -oo  <  x  <  h       (5.25) 

4<k»  €  /    ii       o  J  xnc  —   — 

4   B 


oo 


—         Ih   (X^H  ^[klx-x'ljdx': 

4wf    J   hA     o   l     '     '  J 
-h 


E      for   -h  c   x  <:  oo 
inc  -   — 


and 
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(5.26) 


k 

4W£ 


I   (x^H  ^[klx-x'Mdx 

I    00       o    L         J 
-OO 


Ej     for   -oo  <  x  <  oo 
inc  —   — 


(5027) 


Thus  from  the  above  we  obtain 
h 


40)6 


oo 


f  ih  (x)H  ^^klx-x'Ddx^   I  I,  (xk  (2)(k|x-x'!)  dx' 

J     B     °  J    hA    ° 

-oo  -h 


.oo 


(2) 


I  (x^  H  ^  '(k|x-x'|)dx  a  -EJ       for   -h  <  x  <  h. 
oo     o  inc  —   — 


-oo 


(5028) 


Thus  from  Equation  (5024)  and  (5.28) 


k 

4to£ 


[I    (x')    -   Kx')]   H    ^(klx-x'Ddx  =  A  E, 
so  v 


-h 


x) 


.-h 


k 

4w£ 


(2), 


/.  ,  .    / 


[I,     (x)    -    I    (x^H    ^'[k|x-x'|]dx' 
n  oo  o 

B 


-oo 


oo 

+     j     tXh    (x)    "    I00(^]Ho(2)[k|x-x'!]dx 


(5029) 


for  -h  <  x  <  h„ 


The  above  equation  states   that   the  contributions  due   to   I,     (x;    for  x  '<£   h, 

B         ~ 

Ih  (x)  for  h  <  x'<  oo  and  I  (x^  for  x/  -oo  <   x'<  oo  satisfy  the  boundary 

A  -   -        oo  -   - 

conditions  of  the  problem  exactly.   It  is  interesting  to  note  that  then 
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A  E  i.e.  the  mismatch  in  the  boundary  condition  is  due  to  the  fact  that 

only  parts  of  these  current  1^.  (x  ),    Ih^C*7)  anc*  I  (x  )  in  the  finite 

region  AB  have  been  taken  into  account  for  the  approximate  current 

distribution  I(x).   Thus  it  is  obvious  that  A  E  is  the  contribution  of 

those  parts  of  Iv>»(x  ),  Iv,  (x /)  and  I  (x  )  which  have  been  excluded  for 
A        B  °o 

approximating  I  (x  )  as  given  by  the  Equation  (5.29). 

Since  both  Iv,  (x  )  and  Ih  (x  )  tend  to  I  (x^   for  x   sufficiently 
nA  nB  oo 


length,  both  Iv.  (x7)  -  I  (x')  and  Iy^Cx7)  -  I  (x  )  are  very  small  within 


away  from  the  edge^  when  the  width  of  the  strip  is  comparable  to  a  wave- 

the  range  of  integration  as  given  by  the  right  hand  side  of  the  Equation 
(5.29).  It  is  therefore  expected  that  A  E  will  be  a  very  small  quantity 
especially  when  h  is  very  large. 


P 


%  — H 


FIGURE    7 


As   discussed   before      Iv,      and    !>.      are    of    the    form 


B 


.JT 


2e  ikr   cos  a 

I.     ,_     .  =  !k   sm  a    e 

half   plane 


a 


/2kr   cos  7j 


Uji/ff 


-X     ^  /2k  a    -ikr 

e        dX  -   i  7— sin  -e 
V  7T  2 


i a 

-  /2kr   cos  — 
2 


(5.30) 


l/2kr   cos  — 
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a 


When   the   angle  a   £   TT  and      /2kr   cos   —     >  >  1} 


!  _iX2 
e    d\  can 


/2kr  cos  — 
2 


be  expanded  asymptotically  as  follows:   Let  a  =   2kr  cos  a/2  and 
B  =  /2/7T  A.   Then 


f  -"^Wl  / 


B       .    77      2 

-l   —   t 

e  dt 


B  77    ,2 

r  _i  2 

v/27f  e  dt 


-A 


-B 


Since   Fresnel    integrals   are   odd    functions 


-a2     -F 

|l(i-i)+i§ge    2    -s^^.ori^)] 

7TB  B 


I'  2 


2  a 

—  -i2kr   cos     — 

-  [  (l-i)+e  (- 


2kr   cos 


3  a     / 

4kr   cos      —  /77kr 
2    v 


Therefore 


*  °  <l?>2  1 


.77 
X4 


~  — — 7=1  k   sma,/—   (l-i) 
half    plane  -U(j./7T  If  2 


ikr   cos  a 


+   k   sin  a,/  -  e 


77     -ikr 


i/77kr   cos  —        4kr   cos      —  •    i/77kr 
v  2  2      v 


,  2k       .a      -ikr. 
)    -    l  /  —  sin  —  e 
1  r  2 


77 

i—  .a      -ikr 

4  k   sm  —   e 

2k      .              ikr   cos  a        e  2 

—  sin  a    e  -  — —     — — 

^u  aiu         , /o^-,.„   ,.„ 2   a 


°^         v/27Tkr   kr   cos      - 
V  2 


+  0  (L.)2- 


kr 


(5.31) 

The  current  induced  on  the  surface  of  an  infinite  plane  under  the  same 
condition  is  given  by 


2  ikr  cos  a 

I   =  — -  k  sin  a  e 
oo   Wu 


Therefore 
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.  7T 

1  -        a   -lkr 
4   k  sin  —  e 

j      -  i   =  -  £ 2 

half    °o      u)" 


^    1/2F  (kr)2-  cos^  - 

2 


7T 
1  4 


/27n 


a      .. 
sm  —    -lkr 
2    e 


2  a      3 
COS  2   (kr)2 


where  "H  =  /  tL 


(5.32) 


Here  it  may  be  noted  that  the  phase  is  referred  to  the  edge  of  the 

3/2 
half  plane  and  that  the  magnitude  is  proportional  to  (kr)    . 

Since  r  is  given  by  (h-x')  for  1^  (x7)  and  by  (h+x  )  for  Ih  (x),  within 

B  A 

the  range  of  integration^  r  is  always  greater  than  2h. 

Hence  the  asymptotic  expansion  is  justified  for  Ih  (x  )  if  a  ^  7T 

and   v/4kh  cos  a/2  >  >  1,  and  for  Ih  (x  )  if  a  £   0  and   /4k h  sin  a/2  >  >  1. 

A 


r 


■  H 

1  4   sin  -        -ik(h-x') 
e         2   e 


Ih  (4  "  I  (x')  „   - 

B  -   /SB*     cos2*  [k(h-x^]3/2 


ikh  cos  a 


I   (X7)  -  T   (X')  „  - 


-ikh  cos  a 


Vh 


OC' 


i/2h  sin2  I  [kCh+x7)]3/2 


(5.33) 


Therefore 


AE 


.7T 

k e 

4w£t>,  ~ 


/"   sin  -   -ik(h-x/-h  cos  a) 

— ITa"- 372—  Ho  ^[klx-x^ljdx 

J.   cos2|   [k(h-x')]3/2 


00      a  / 

f  cos  —   -ik(h+x  +h  cos  a) 

I     —A „,»      H(2)[k|x-x  |]dx' 


2  a 


h    sin  -   [k(h+x/)] 


3/2 


46 


AE  = 


.k  e 


4ooc  n 


«   277 


/•      sin  -        -ik(h+x'-h   cos  a) 

—£KZ_ Ho(2)[k(M/)]dx/ 

h      C°S     2  [kCh+x^]372 


f     cos  —        -ikCh+x^h  cos  a)         >„* 

+  J   TIT*6  rkfh   ,..3/2     ■,'w'-"" 

<       sin     -  [k(h+x')] 


(5.34) 


Using  the  integral  representation  of  Hankel  Function  given  by 


OO 


H     (x) 
o 


-ixt 


{  JT^i 


dt 


we  obtain  from  Equation  (5.34) 

7T 
„2  X   4 

4u£*l      y=r     TTx 


a  °°       °°  / 

f         f       -ik(h+x  -h   cos  a)    -    ik(x+x   )t 
e 


^x  r  r 

2  a 


cos 


\   {        [Mh+x')]372  /7-i 


dt   dx 


cos  — 


a  oo       oo 


+   —2^         I  J 

sin     —      J        J 

2        h        1 


f        r       -ik(h+x4-h   cos  a)    -   ik(x'-x)t 


[k(h+x')]3/2  /j 


t2-x 


Now  consider 


dt   dx 


(5.35) 


oo        oo 


oo  /  / 

t        -ik(h+x  )    -    ik(x  +x)t 

h 


{    [k(h+x')]3/2  7 


t  -1 


oo  oo  / 

f        +ikxt  r        -ik(l+t)x 

-ikh  e—  ,.  /       e  / 

_      e  _______  dt  — — — — — -—-  dx 


>x     \ft   -1 


Jh      [k^+x7)]372 


by   changing   the   order  of   integration, 
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But 

°e-ik(i+t)x/    ^       2    ^_ikh(1+t)    «  ^  r     r  (r+i) 


dx  «i    -  ~     e  2  (i) 


J.       [kCh+x')]372  k  r=l  (l+t)r(2kh)r+^ 


by  integration  by  parts  when  kh  >>1. 

Pe-ik(h+x/)-ik(X±X)t   HtH    •  2      -2ikh  *  n    r  F  (r+j)         fV^h+x) t 

dtdx^,-  t-  e  S  (i)      — — T — — —  dt 


JhJ       [k(h+x')]3/2   i/t2^  k  r=l  (2kh)r+i    J     (l+t)r(/?^i 

(5.36) 
But 

r°°     -ik(h+x)t  /°     -ik(h+x)z 

I        e  =      e-lk(h+x)         I      e  dz 

1       (l+t)r  /t^l  J        (2+z)n+V^ 

I  o 

by  change  of  variable. 

When  the  point  of  observation  is  such  that  k(h  +  x)  is  sufficiently 

large  to  permit  asymptotic  expansion, 

r°     -ik(h+x)z 
£_ Z^dz 

J   (2+z)n+Vz~ 
o 

can  be  evaluated  asymptotically.   If  h  +  x  =  0,  the  above  can  be  evaluated 

exactly,  while  in  general,  if  k (x  +  h)  is  not  large  enough  to  permit 

asymptotic  expansion,  the  upper  bound  of  the  above  can  be  determined. 

Thus  the  following  three  cases  will  be  considered. 

Case  I.   The  point  of  observation  is  at  or  near  the  center  of  the  strip 

so  that  k(h  +  x)  is  sufficiently  large. 

Consider 

f   -ik(h+x)z 

",   dz 
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-Z. 


v> 


Z  -Tla-ne 


^U 


FIGURE  8 


taken  around  the  contour  C  in  the  complex  z  plane  as  shown  in  Figure  8^ 
in  which  the  branch  points  are  at  -2  and  0,   Let  a  branch  cut  be  made 

along  the  real  axis  from  u  =  0  to  u  =  «  and  to  u  =  -2  and  let  the  branch 

1 
of  \fz   be  chosen  in  which  the  value  of  fz   is  i | z i  2  when  z  is  real  and 

negative.   The  integrand  is  analytic  along  the  contour  C  as  shown  in 

Figure  8  and  hence  the  value  of  the  integral  around  the  path  is  zero. 

The  integral  around  the  large  circle  tends  to  zero  when  the  radius  of 

the  arc  tends  to  infinity.   The  integral  around  the  small  circle  also 

goes  to  zero  as  the  radius  of  the  small  one  tends  to  zero. 

Hence 


r       -ik(h+x)z        -i-r 
e     —     j        4 
x    dz  =  e 


OO 


(2+z)r- 


\  z 


-k(h+x)v 
e    — 

(2-iv)  +2  {v 


dv 
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which  by  Watson's  lemma  can  be  expanded  asymptotically  and  is  given  by 


OO 


-ik(h+x)z 
e     — 

(2+z)r+i  k 


H 


_i—  co  F  (m+J)a 


dz 


r.m 


m=0  [k(h+x)] 


T  ) 


with  h  +  x  >  0 


(5.37) 


where  a   is  defined  by 

rm 


S  a   v"  =  (2-iv)-(r+i) 
rm 
nt=0 


(5.38) 


To  sum  up  the  result,  for  x  «  0;  such  that  k(h  +  x)  ^>1, 


<AE  ~  - 


-i3kh 


T!   ./27T 


-ik(x-h  cos  a)  £    r  T  (r+|)   *   r(m+*)arm 

e  2  (l)   — - a— T  2 


cos 


2  a 


r=l     (2kh)r+2  m=0  [k(h+x)]m+2 


a 


OO 


„  r,   u    °o   r  (m+J)a 
cos  2   lk^x-h  cos  a)  ^,  , .  .r     r(r+.,)    „  rm 

+  -air  e  S  (i)     -^r-  _., 


sin 


r=l      (2kh)   2  m^O   [k(h-x)] 


It  may  be  noted  that  since  in  the  expression  for  current  along  a 

_5 
half  plane  as  given  by  Equation  (5.31)  is  accurate  to  the  order  of  (kr)  2? 

the  above  equation, (5.39)  will  yield  results  accurate  up  to  the  same 

order.   However  by  considering  terms  involving  higher  powers  of  (— ) 

'  kr 

in  Equation  (5.31),  Equation  (5.39)  may  be  made  to  yield  results  accurate 

up  to  any  higher  order  of  (7—) . 

kr' 

Case  II:   The  point  of  observation  is  at  the  edge  i.e.  x  -  +  h. 
In  this  case 


r  -ik(h+x)t 

e     — 


r  r3>" 
(1+t)   /t  -1 


=  dt 
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r(r+-|)  „  *      rm 

(2kh)r+*  m-0    (2kh)m+* 


r-1) 


2kh) 


r7* 


(5„40) 


P(1f  ^JiziU^ 

1  2r'2kh)r+2 


°o     r(m+|)a 


rm 


^2  HteO    (2kh)m+2  _ 


(5.41) 
edges   say  B  ,    i.e. 
ge   enough   to   justify 


O    PLEASE  PEEL  OFF  AT  EXTENDED  SHEETS  AND  SEND  COPIES 

^    1   AND  2  (HELD  TOGETHER)  TO  BINDERY  NUMERICALLY  ARRANGED. 

in  terms  of  k(h-x).  Then  the  result  obtained  in  Case  1  for  the  second 
terra  in  the  right  hand  side  of  Equation  (5,39 »  is  not  valid.  However, 
in  this  case  an  upper  bound  of  AE  can  be  obtained  as  follows; 


cc 


f        -ikfh-xit 
e 

/t^l  (1+t)1 


dt 


2r   r(r+A) 


(5.42) 
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Therefore 


AE 


at   x=h 


Tl    /27T 


sin 


2        -ikh(4-cos  a) 


cos 


2  a 


2  CD 
r=l 


r   r(r+i) 


co     T(m+h)a 


rm 


(2kh) 


r+2   m=0    (2kh) 


m+, 


while 


AE 


(r-1)! 


.IT  a 

1 —    COS    —  oo 

+   e  4  2  e-ikh(2+cos  a)      z  (i)r 

sin2  \  r-1  2r(2kh)r+" 


at   x--h 


Tl    /27T 


a 

sin   2 

2  a 

cos     — 

2 


3T-. 
-i      kh(2-cos  a)-  -jj-  oo 

:  ''  S  (i) 

r   1 


15.40) 


(r-1)  I 


2    (2kh) 


•,w„  ->      °°  r/       ,v  °o     r(m+A)a 

s-ikh(4+cos  a)     s  (i)r  r(r+&  £  rm 

r=l  C2kh)r^"2   BteO    (2kh)m+2 

(5.41' 

Case  III-   The  observation  point  is  near  one  of  the  edges  say  B,  i.e. 
h  and  x  are  of  the  same  order  and  k(h-x)  is  not  large  enough  to  justify 
asymptotic  expansion  of 


r°°     -k(h-x)v 

e 


dv 


(2-iv)  +2  {v 
o 

in  terms  of  k(h-x).   Then  the  result  obtained  in  Case  1  for  the  second 

term  in  the  right  hand  side  of  Equation  (5,39)  is  not  valid.   However, 

in  this  case  an  upper  bound  of  AE  can  be  obtained  as  follows: 


C        -ik(h-x'»t 


/t  -1  (1+t) 


dt 


i_   C>-l)j 
2        n'r+A) 


(5,42) 
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Now 


00         oo 


^-kl  (h+x  ')-ik(x'-x)  t 

k(h+x')]3/2  /t2"l 


oo 


dt    dx' 


2      -2ikh  J*3,  .r   T  (r+i) 
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2   e 
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4  -ik   <h(3+cos   a)-x>     - — ^ —  <    2j    — - —         ^  -~   +   i     2, 
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x2r+3  f 

>  J 


(5.43) 


Similar   results   can   be   obtained   for   the  observation  point   x  near   the 
other   edge  A,    i.e.,    x  «  -  h.      In   that   case 
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It  as  interesting  to  note  that  the  first  term  in  all  these  asymptotic 

1   2 
series  involves  at  lease  (  --)   and  succeeding  terms  are  of  higher  order 

kh 

in  (s-r).   Thus  if  kh  is  sufficiently  large,  a  good  approximation  can  be 

obtained  even  if  only  the  first  few  terms  of  the  series  are  considered, 

A  E 
In  a  later  section,  A  E  expressed  ab  — X  100  for  various  h  and 

E  . 
mc 

a  have  been  computed  and  plotted  in  order  to  determine  the  extent  of 
deviation  from  the  boundary  condition  with  assumed  current  distribution. 
It  is  found  that  this  deviation  is  very  small  for  the  sizes  of  strips 
used  in  the  construction  of  the  zoned  mirror  in  which  the  width  of  each 
zone  is  of  the  order  of  a  wavelength  or  more,   Thus  the  use  of  the 
approximate  current  distribution  for  the  computation  of  the  field  of 
the  zoned  mirror  assembly  has  been  justified. 
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6    NUMERICAL  COMPUTATIONS 

For  the  purpose  of  computation  the  dimensions  chosen  are: 

(a)  F  -  f/D  -  ,556    where  F  Is  the  focal  length  and  D  is  the  diameter 

(b)  f       =  10\.    of  the  aperture. 

(c)  Total  number  of  Zones  =  11. 

Equations  (5,8),  (5,9)  and  (5.10)  are  used  for  computation  of  I . .  while 
Equation  (5,16)  is  used  for  computation  of  coupled  current  I. ..   The 
expression  for  I  used  for  the  above  are  given  by  Equations  (5.13),  (5,14) 
and  (5.15).   Equations  (5,15,  (5.17)  and  (5,18)  are  not  used  because  when 
the  arguments  of  the  Bessel  functions  of  the  two  kinds  approach  each  other, 
the  convergence  of  the  series  given  by  Equation  (5.19)  becomes  very  slow 
and  hence  the  number  of  terms  of  the  series  required  for  a  certain  accuracy 
becomes  extremely  large. 

Once  the  current  distribution  across  a  given  strip  has  been  calculated 
from  the  above,  the  electric  field  at  any  point  due  to  each  strip  is  com- 
puted from  Equation  (3.19),   In  computing  the  field  due  to  any  strip,  the 
asymptotic  form  of  Hankel  function  is  used.   Since  the  geometry  is  such 
that  the  argument  of  the  Hankel  function  in  Equation  (3,19)  i.e  ,  kr  is 
much  greater  than  unity. 

The  first  zone  is  a  smoorh  parabola  of  finite  dimension.   The  rigorous 

solution  for  this  is  unknown.,   However,  if  the  parabolic  zone  is  replaced 

by  a  flat  strip  of  the  same  aperture,  it  is  possible  to  compute  I. .,  I   and 

fields  due  to  I    and  1    +  1   .   It  is  found  that  not  only  was  I    very 
ii       ii     ij  ij 

small  compared  to  I   ,  but  the  fields  due  to  I    and  I   +.  I .  .  were  practi- 

ii  ii       ii     ij 

cally  the  same.   Further,  we  can  reasonably  assume  that  I. .  of  the  first 

zone  either  as  a  flat  strip  or  a  smooth  parabola  should  be  of  the  same 

order  of  magnitude.   Therefore  we  conclude  that  i  .+1  .  and  I  .  produce  almos 

n   ij       n 

the  same  field  for  a  smooth  parabola  also.   However,  I  .,  the  current  induced! 

'       ii 
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on  a  finite  smooth  parabola  due  to  an  incident  plane  wave  is  still  unknown; 
but  it  is  generally  known  that  the  physical  optics  solution  for  a  finite 
parabola  agrees  very  closely  with  experimental  results.   Further  according 


Table  1 

Comparison  of  the  contribution  of  the  first  zone 

(flat  strip)  to  the 

image 

pattern  (a)   with  physical  optics  solution    (b) 

with  edge  effect. 

Magnitude  of  the  E  Field 

y 

Case  (a) 

Case  (b) 

0 

1,29 

1,28 

-.2 

1„27 

1.26 

-.4 

1.22 

1,21 

-.6 

1.14 

1.14 

-.8 

1.05 

1.05 

-.9 

1.01 

1.01 

-1.0 

,97 

.96 

to  the  computation  (Table  1)   the  physical  optics  solution  for  the 
corresponding  fiat  strip  does  not  differ  significantly  from  that  including 
the  edge  effect,  since  the  strip  is  quite  wide  in  terms  of  wavelengths 
(6.245  A-) .   Therefore  it  seems  reasonable  to  use  the  physical  optics  solution 
for  the  smooth  parabola  representing  the  central  zone.   In  other  words, 
the  significant  contribution  of  the  first  strip  is  that  due  to  the  geometri- 
cal optics  current  regardless  of  whether  the  first  zone  is  a  flat  strip  or 
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a  portion  of  a  smooth  parabola,   However,  this  does  not  imply  that  the 
contribution  due  to  the  flat  strip  is  the  same  as  that  due  to  an  equivalent 
smooth  parabola.   This  becomes  clear  if  one  observes  that  in  the  present 
case  the  geometrical  path  difference  at  the  edge  between  the  two  is  of  the 
order  of  \'4<       The  difference  in  over-all  contribution  is  indicated  in 
Figure  9  for  normal  incidence,   It  shows  that  in  actual  construction  the 
first  zone  should  not  be  approximated  by  a  flat  strip. 
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2-0 


FIGURE  9   CONTRIBUTION  OF  IMAGE  PATTERN  ON  THE  FOCAL  PLANE  OF 
THE  ZONED  MIRROR  DUE  TO  THE  CENTRAL  ZONE  ONLY  FOR 
NORMAL  INCIDENCE 
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7.   DISCUSSION  OF  RESULTS 

7.1   Error  in  the  Satisfaction  of  the  Boundary  Condition  Due  to  the 
Approximation  of  the  Current  Distribution  on  the  Strip 

The  magnitude  of  the  error  in  the  satisfaction  of  the  boundary 

condition.!  A  E  I ,  on  the  surface  of  an  isolated  strip  due  to  the  approximate 

current  distribution,  expressed  in  percentage  of  incident  electric  field, 

is  plotted  as  a  function  of  the  observation  point  on  the  strip  in  Figure  10 

for  strip  widths  of  2\}     X,  .75  X  and  0    .  5X  with  normal  incidence  and  in 

Figure  11,  I  A  E  I  is  plotted  for  a  strip  of  width  X   for  incident  angles  of 

o     o     o        o  II 

0  ,  15  ,  25   and  45  .   The  graph  shows  that  I  A  E  I  varies  from  point  to 

point,  across  the  strip.   The  average  of  I  A  E  I  ,  however,  increases  if  the 
strip  width  decreases  or  for  a  given  width  if  the  angle  of  incidence 
increases.   For  the  dimensions  of  the  strips  used  in  the  zoned  mirror  the 
average  I  A  E  !  is  less  than  about  0,5%  for  normal  incidence;  and  about 
2.0%  for  oblique  incidence.   It  is  therefore  clear  that  the  half  plane 
current  distribution  is  a  very  close  approximation  of  the  exact  distri- 
bution,,  It  is  doubtful  whether  much  greater  accuracy  could  be  achieved 
by  using  the  exact  solution,  due  to  the  difficulty  of  computation  as 
discussed  in  great  detail  in  the  appendix. 
7„2  Current  Distribution  (1,^  and  Ij  •) 

The  surface  current,  density  I.,  on  an  isolated  strip  due  to  an 

incident  plane  wave  and  the  first  order  surface  current  density  I.  . 

•^  J 

induced  on  the  strip  due  to  the  presence  of  other  strips  have  been 
plotted  for  different  configurations  of  the  strips  and  different  angles 
of  incidence  in  Figures  12  to  17.   The  curves  of  I. .  show  the  correct 

IX 

singularity  at  the  two  edges  as  expected  since  they  are  ha If -plane 
solutions.   The  current  drops  sharply  to  a  value  well  below  that  at 


_o —  1=.5X 
—a—  Jt-.75X 
— □ —  J-X 


58 


FIGURE  10   I  A  E I  IN  PERCENT  OF  INCIDENT  FIELD  VS  POINT  OF 
OBSERVATION  ANGLE  OF  INCIDENCE  a  =■  0 
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Angle  of  Incidence 


O  c:  -  o 

A  cx^  15° 

O  CX  =25° 

O  ck-   45° 


FIGURE  11   lA  E I  IN  PERCENT  OF  INCIDENT  FIELD  VS  POINT  OF 


OBSERVATION  (STRIP  WIDTH  =  IX) 
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FIGURE  12    SELF  CURRENT'  AND  l MUTUAL'  CURRENT  ON  STRIP  6 
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FIGURE  13   SELF  AND  MUTUAL  CURRENT  ON  STRIP  5 
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FIGURE  15   »  SELF'  AND  ■  MUTUAL'  CURRENT  ON  STRIP  3 
for  a  ~:  0  and  15 


64 


Up 


FIGURE  16   '  SELF'  AND  '  MUTUAL'  CURRENT  ON  STRIP  2 
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FIGURE  17   SELF  AND  MUTUAL  CURRENT  ON  STRIP  1 


66 

the  middle  as  the  point  of  observation  is  moved  away  from  the  edge.   It 

then  increases  again  approaching  the  geometrical  optics  value  near  the 

middle.   The  rate  of  decrease  is  reduced  as  the  angle  of  incidence  is 

increased. 

The  first  order  coupled  current  distribution  I.   also  has  similar 

ij 

singularities  at  the  two  edges  as  evident  from  Equation  (5.13).   This 

current  distribution  varies  with  the  size  of  the  adjacent  strips,  their 

relative  inclination  and  their  distance  from  the  strip  under  consideration,, 

In  general  I .  .is  in  the  order  of  10%  of  I   .   Hence  the  second  order 
ij  ii 

coupled  current  will  be  in  the  order  of  1%  of  I...   Thus,  although  it  is 
possible  to  obtain  the  second  order  solution  in  the  present  iteration 
procedure,  it  seems  to  be  insignificant  and  it  is  probably  within  the 
error  introduced  by  substituting  the  exact  solution  of  an  isolated  strip 
by  half  plane  solution  in  actual  computation  as  explained  before. 
7.3   Image  Patterns 

To  demonstrate  the  difference  in  performance  with  various  orders  of 
approximation,  the  image  pattern  of  the  zoned  mirror  has  been  studied  with 
the  following  current  distributions  on  the  strips  due  to  a  plane  wave  at 
various  incident  angles: 

(1)  The  geometrical  optics  current  distribution  i,e.  the  current 
distribution  on  the  strip  which  is  obtained  as  if  it  were  a  part  of  an  infinite 
plane.   In  this  case  the  effects  of  both  the  edges  and  coupling  between 
strips  are  neglected. 

(2)  A  current  distribution  which  takes  into  account  the  edge  effect 
but  ignores  the  effect  of  coupling  between  strips  i.e.  the  zero  order 
solution  of  the  integral  equation  given  by  Equation  (3.15). 
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(3)   A  current  distribution  which  considers  both  the  edge  effect  and 
the  first  order  coupling  between  adjacent  strips  as  given  by  Equation  (3.17) 

All  computations  are  done  by  ILLIAC,  the  electronic  high  speed  computer 
at  the  University  of  Illinois.   The  time  required  for  computing  an  image 
pattern  with  one  angle  of  incidence  for  three  different  current  approxi- 
mations is  roughly  given  below; 

(a)  Case  1     14  minutes 

(b)  Case  2     16  minutes 

(c)  Case  3      1  hour 

For  comparative  study,  the  image  pattern  of  an  equivalent  smooth  parabola 
i.e.  a  smooth  barabola  with  the  same  focal  length  and  aperture  area  as  the 
zoned  mirror,  is  also  computed,,   These  are  all  shown  in  Figures  18  to  23, 

With  normal  incidence,  the  maximum  intensity  due  to  a  smooth  parabola 
and  Case  1  of  the  above,  are  almost  identical,  while  Cases  2  and  3  yield 
almost  the  same  value.,  which  is  about  four  percent  less  than  that  of  a 
>oth  parabola.   Thus  the  edge  effect  reduces  the  gain  slightly. 

The  secondary  maxima  occur  at  the  same  positions  for  the  smooth 
parabola  as  for  the  Case  1  and  Case  3  of  the  zoned  mirror  while  for  Case  2 

are  slightly  displaced.   The  value  of  the  secondary  maxima  is  the 
larges       the  smooth  parabo]      Le  for  Cases  1,  2  and  3  of  the  zoned 
mirror  they  are  the  same,  having  a       approximately  60  percent  of  that 
the  smooth  parabola. 

As  the  angle  of  incidence  is  changed  from  normal  the  maximum  intensity 
falls  off  rapidly  for  the  smooth  parabola.   Case  1  also  shows  a  drop  while 
for  Cases  2  and  3  the  drop  is  very  small  until  the  angle  of  incidence 
exceeds  15  .   The  locations     '  e  primary  and  secondary  maxima  for  all 
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FIGURE  22   IMAGE  PATTERN  ON  THE  FOCAL  PLANE  a    =  20 
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the  three  cases  of  the  zoned  mirror  remain,  the  same  while  those  of  the 
smooth  parabola  change  considerably.   The  sharpness  of  the  intensity 
distribution  curves  (i.e.,  the  focussing  ability)  for  the  smooth  parabola 
drops  off  considerably. 

The  ratio  of  the  magnitude  of  the  two  secondary  maxima  for  the  three 
cases  of  the  zoned  mirror  remain  constant  and  equal  to  unity  up  to  an  inci- 
dent angle  of  about  20   while  it  changes  rapidly  for  a  smooth  parabola 

as  the  angle  of  incidence  is  changed,  e.g.,  a  change  in  the  angle  of  inci- 

o 
dence  from  0  to  15  changes  the  ratio  from  1  to  4.25.   (See  later  in 

Figure  26) . 

The  absolute  value  of  the  secondary  maximum  for  the  zoned  mirror  shows 
a  slight  increase  as  the  angle  of  incidence  is  changed  from  normal.   The 
secondary  maxima  for  Case  1  deviate  more  and  more  from  those  of  Cases  2 
and  3  as  the  angle  of  incidence  is  increased  and  with  an  incident  angle 
of  20  .  the  two  secondary  maxima  for  Case  1  give  a  ratio  of  1.2  as  against 
1.01  for  Case  2. 

The  fact  that  the  intensity  distributions  of  the  electric  field  for 
Cases  2  and  3  are  very  nearly  the  same,  indicates  that  the  contribution 
of  the  coupled  current  to  the  field  is  negligible. 

These  results  indicate  an  important  fact  that  the  coupling  effect 
between  strips  can  be  ignored  without  sacrificing  the  accuracy.   It  may 
also  be  remarked  that  the  coupling  effect  which  is  the  result  of  the 
first  order  iteration,  is  the  most  time  consuming  step  in  the  computation 
as  indicated  in  the  beginning  of  this  section. 

Though  there  is  some  difference  between  Case  1  and  Case  2  particularly 
for  oblique  incidence,  the  physical  optics  solution  as  given  by  Case  1 
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yields  a  reasonably  good  approximation  (though  at  less  accuracy)  even  with 
22  discontinuities  due  to  the  edges  in  the  present  system.   The  simpli- 
fication of  computation  is,  however,  considerable. 

It  is  interesting  to  note  that  the  two  secondary  maxima  on  either 
side  of  the  principal  are  more  symmetrical  for  Cases  2  and  3  than  for 
Case  1  indicating  that  the  effect  of  the  edges  seems  to  reduce  the  "coma1 . 

According  to  parageometrical  optics  the  focal  plane  is  actually  not 
the  best  image  surface.  Accordingly  an  image  pattern  is  obtained,  as 
shown  in  Figure  24,  in  a  plane  (for  simplicity  of  computation  only  a  plane 
surface  is  considered)  passing  through  the  best  image  point  and  perpendicular 
to  the  focal  axis,  for  15  incident  angle.  It  show  that  the  gain  is  slightly 
higher  but  the  minor  lobes  become  large  and  more  asymmetric. 
7.4  Radiation  Patterns 

The  curves  in  Figures  25,  26,  27,  28,  29  and  30,  which  are  derived 
with  the  aid  of  the  reciprocity  theorem  from  the  intensity  distribution 
curves  discussed  before,  show  the  important  and  interesting  properties  of 
the  zoned  mirror  as  an  antenna  in  comparison  with  a  smooth  parabola,  with 
isotropic  feed   It  is  clear  from  Figure  25,  that  for  a  zoned  mirror  the 
deflection  of  the  beam  of  the  radiation  pattern  from  the  axis  is  propor- 
tional to  the  displacement  of  the  feed,  while  it  is  not  so  for  the  smooth 
parabola.   For  the  same  scan  angle,  a  larger  displacement  of  the  feed  is 
required  for  the  smooth  parabola  than  for  a  zoned  mirror  and  this  difference 
in  displacement  increases  as  the  scan  angle  is  increased.   Figure  26  shows 
the  maximum  field  intensity  on  the  focal  plane  for  various  angles  of  inci- 
dence while  Figure  27  shows  the  variation  of  the  directive  gain  against 
the  scan  angle.   The  gain  of  the  zoned  mirror  is  substantially  constant. 
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FIGURE  27   RELATIVE  GAIN  OF  THE  RADIATION  PATTERN  VS  SCAN  ANGLE 
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FIGURE  28   BEAM  WIDTH  OF  THE  RADIATION  PATTERN  VS  SCAN  ANGLE 
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FIGURE    30      MINOR   LOBE   LEVEL  AT   THE   RADIATION   PATTERN  VS    SCAN 
ANGLE 
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For  a  15   scan  angle,  the  loss  in  gain  is  only  0.27  db  as  referred  to  the 
on  focus  feed  case,  while  for  a  parabola  the  loss  is  1.36  db„ 

Figure  28  shows  that  the  change  in  beam  width  for  0   to  15   scan  angle 
is  roughly  30%  for  the  smooth  parabola  against  5%  for  the  zoned  mirror. 
Figure  29  is  a  plot  of  the  magnitude  of  the  first  two  secondary  maxima  of  the 
image  pattern  against  their  positions  in  the  focal  plane.   In  Figure  30, 
a  plot  of  minor  lobe  level  of  the  radiation  pattern  against  the  scan  angle, 
it  may  be  noted  that  the  minor  lobe  level  is  smaller  for  the  smooth 
parabola  than  that,  of  the  zoned  mirror  for  the  on  focus  feed  case.,  the 
figures  being  -16,14  db  and  -14.43  db  respectively.   This  is  expected 
since  the  zoned  mirror  having  reflecting  strips  at  a  constant  distance 
from  the  focus,  is  uniformly  illuminated  while  the  smooth  parabola  is  not. 
The  minor  lobe  level  of  the  radiation  pattern  goes  up  rapidly  as  the  off- 
setting of  the  feed  is  increased  for  the  smooth  parabola  whereas  it 
increases  very  slowly  for  the  zoned  mirror.   For  example,  as  the  scan 
angle  varies  from  0   to  15   the  minor  lobe  level  is  raised  by  9.03  db  for 
the  smooth  parabola  and  only  1,95  db  for  the  zoned  mirror.   At  15   scan 
angle  the  minor  lobe  level  is  only  -7.2  db  for  the  smooth  parabola  and 
-12.6  db  for  the  zoned  mirror. 

A  more  interesting  result  so  far  as  the  coma  correction  is  concerned, 
is  shown  in  Figure  31,   It  is  seen  therefrom  that  the  ratio  of  the  first 
two  secondary  maxima,  (a  function  of  coma  aberration)  is  substantially 
constant  for  a  zoned  mirror  in  contrast  with  that  for  a  smooth  parabola. 
It  clearly  indicates  that,  the  zoned  mirror  like  the  one  considered  is  very 
effective  in  the  reduction  of  the  coma  aberration. 
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8.   CONCLUSIONS 

8.1  The  Approximate  Current  Distribution 

A  study  of  the  I  A  eI  curves  reveals  that  the  approximate  solution  is 
very  close  to  the  exact  one.   Because  of  the  ease  of  computation,  it  has 
more  practical  value  than  the  exact  solution.   The  computation  for  the 
exact  solution  may  be  extremely  difficult  if  not  impossible. 

8 . 2  The  Technique  of  Analysis 

In  spite  of  a  large  number  of  discontinuities  in  the  structure,  the 
generally  used  physical  optics  solution  offers  a  sufficiently  accurate 
result.   This  method  is  simple  and  straightforward.   For  higher  accuracy 
the  edge  effect  should  be  considered  but  the  coupling  effect  seems  to  be 
of  less  importance 

In  view  of  the  fact  that  the  contribution  of  the  edge  effect  is  only 
a  higher  order  correction  to  the  physical  optics  solution,  it  is  expected 
that  the  computation  of  this  contribution  may  be  obtained  simply  by 
replacing  each  edge  by  a  line  source  of  proper  intensity  in  which  case 
the  calculations  involving  Fresnel  integrals  can  be  avoided. 

8.3  The  Performance  of  the  Zoned  Mirror 

The  effectiveness  of  a  coma-correlated  zoned  mirror  has  been  the 
object  of  controversy  for  some  time.   The  present  investigation,  made 
on  a  two-dimensional  cylindrical  zoned  mirror  which  is  free  from  astig- 
matism, enables  us  to  evaluate  its  effectiveness  in  correction  of   coma" 
particular] y  for  a  system  with  a  small  value  of  f/D  (the  smallest  F-number 
of  this  system  is  1/2),  and  uniform  illumination,  without  complications 
due  to  other  aberrations. 
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A  comparison  with  the  performance  of  an  equivalent  smooth  parabola 
has  indicated  that  the  zoned  mirror  with  small  F-number  and  unifrom 
illumination  has  a  performance  much  superior  to  that  of  a  smooth  parabola 
in  the  reduction  of  the  'coma'1  aberration.   The  beam  width,  gain  and 
minor  lobe  level  of  the  radiation  pattern  and  also  the  ratio  of  the  two 
first  minor  lobe  levels  of   the  image  distribution  in  the  focal  plane 
(a  function  of   coma   aber;ation)  are  relativel}'  constant  with  a  zoned 
mirror  for  a  scan  angle  Up  to  25  ?  while  for  a  smooth  parabola  the  ratio 
of  the  first  minor  lobe  level  increases  drastically.   Further  coma 
correction  seems  to  be  possible  by  positioning  the  zones  slightly  off 
the  coma  circle  such  that  the  third  order  and  the  fifth  order  coma 
partially  cancel  each  other. 
8.4   Suggestions  for  Future  Investigation 

A  detailed  numerical  study  of  the  performance  of  zoned  mirror  with 
f/D  ~  1/2  but  with  different  values  of  f/^  followed  by  experimental 
investigation  is  suggested.   Since  the  zoned  mirror  is  a  frequency  sensi- 
tive device,  it  will  be  interesting  to  evaluate  its  chromatic  aberration 
by  studying  the  reduction  in  gain  and  variation  in  minor  lobes  as  the 
operating  frequency  changes. 
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APPENDIX  A 

DISCUSSION  ON  THE  DIFFICULTIES  OF  EXACT  SOLUTION  FOR  NUMERICAL  COMPUTATION 
FOR  WIDE  STRIPS 

It  has  already  been  indicated  that  the  exact  solution  of  the  above 
problem  is  obtained  in  ter  is  of  Fourier  series  expansion  of  various 
Mathieu  functions  as  given  by  Equations  (5.1). 

(a)   The  Fourier  expansion  of  se(s,x)  i?e 


oo 
se  (s,x)  =  L  Deni    cos  (2k+p)x,    where  p  =  0  if  n  is  even       (A„l) 
n        t-n   2k+P  ' 

p  -■    1  if  n  is  odd 


is  no  longer  useful  for  computation  when  s  is  sufficiently  large.   It  can 
be  shown  from  the  asymptotic  behavior  of  the  coefficients  De„,   ,  that 

they  are  extremeljr  large^  even  though  Se  (s,x)  is  small.   For  example 

24 

with  large  s,  we  have7  asymptotically^ 


De  (s)  ~  e  "/s/2s1/4    ■/ ~j     and  I  De  (s)/  I  De  (s)l  «  2  (A.  2) 

based  upon  a  normalization  which  requires  that 

oo 

£  De     =  1 

k=o   2k+p 

For  s  -   2500  corresponding  to  h  approximately  equal  to  8^-  it  would 
be  necessary  to  know  the  coefficients  up  to  22  significant  figures  before 
the  normalization  condition  is  satisfied  up  to  two  significant  figures, 
The  same  difficulty  will  arise  for  functions  of  higher  orders.   This 
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difficulty  can  not  be  avoided  by  using  a  different  normalization.   For 
the  ratio  Se  (s,ff/2)  /  Se  (s,0)  is  independent  of  the  normalization  and 
unless  the  coefficients  were  computed  with  extreme  accuracy,  the  order 

of  the  magnitude  of  the  function  would  remain  undetermined  over  a  large 

8 
part  of  the  S-x  plane.   For  example,  Goldstein  obtained  22  coefficients 

for  S  =  6400  for  five  decimal  places  for  computing  the  function  ce(s,x) 

12 
which  is  proportional  to  Se  (s,x).   But  although  ce(s,JT/2)  can  be  readily 

determined  from  his  coefficients,  the  latter  yield  merely  ce(s,x)  =  0 

for  x  \  55  .   Thus  with  trigonometric  series,  the  individual  terms  are 

too  large  and  too  slowly  convergent  to  be  of  any  use  for  computational 

purposes  when  s  is  large.   The  asymptotic  expansion  seems  to  be  the  only 

feasible  method  for  practical  purposes  under  this  condition. 

The  asymptotic  expansion  of  the  function  Se  (s,x)  is  given  by 

Sen(s,x)  ~  (7rtW1{Po<:x)  "  Pi(x)j  +  W2  (P0(X)  +  Pi(x)f  1  for  "  n/2^   X  v  */2 

(A.  3) 


where 


/s'sin  x    r  ,n        7T  ,2n+l 

Wn    =  e  [cos    (-  +  -)] 


(A.  4) 


.  n+1 
(cos   X) 


-  s/s    sin   xr  .n        7TV  ,2n+l 

W?   =   e  [cos    (--  +  -)] 

,  .n+1 

(cos   x) 


(A.  5) 


_i 
4    /s  32s  256    s3/z 
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_    ,   .  .  v  i       v4    +   86v2    +    105         vA    +   22V2    +57 

PQ(x)  =  1  +  _ ~  +  -  [ ]  +  ...(A. 7) 

4  Vs  cos         512  cos   x        512  cos   x 


(v3  +  3)  sin  x       Sin  X       3         4^   +  44^ 
P^x)  «V +  JJ  S*n  X   ■ g—  fvJ  +  3v  +  ±__+_*i-]  +...  and  v  =  2n+l 

16  V s  cos   x    128s  cos   x  cos   x 

(A, 8) 

The  above  asynptotic  expansion  gives  best  results  when  x  is  very 

small.   The  accuracy  falls  off  as  x  increases.   It  can  be  shown  that  for 

reliable  results 


2  V  v 
COS  X  >  — TTT'  (A. 9) 

—   1/4 
s 


Thus  for  a  given  s  and  x  the  order  up  to  which  the  function  can  be 
evaluated  asynptotically  with  reasonable  accuracy  is  fixed.   For  example^ 
with  n  -  0  i.e  x  ~    1  and  s  =  2500,  Se  (s,x)  can  be  obtained  satisfactorily 
up  to  a  value  of  x  given  by 

2  o   / 

cos  x  >  —  i ,  e  x  -  73  30 


v/50 


For  larger  values  of  x  the  results  obtained  from  the  asymptotic  series 
given  by  Equation  (A3")  will  not  be  satisfactory. 

The  expression  given  by  Equation  (A. 3)  fails  altogether  for  x  =  90  , 

24 
The  above  expression,  however,  can  be  extended  in  the  following  way    to  be 

useful  when  x  -  90  . 

IT 

Let  x  --  -  -  P  be  that  value  of  n  where  the  above  asynptotic  expansion 

IT 

gives  accurate  results.   Now  consider  the  Taylor  expansion  around  x  =  — . 

Let  y  (s.x)  represent  Se   (s,x)  and  y  (s,x)  represent  Se     (s,x). 
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Hence   y^    (s,2T/2)    =   y2(s,JI/2)    =  0   and  y   (s,x),    y    (s,^are   even   and  odd 
functions   respectively  of   V/2   -   x. 
Let 

y^s,    1/2    -  p)    =   yx(s,    IT/2)    T^s^) 

and 

y2(s,    7T/2   -  p)    =   ~y^(s,ir/2)    T2<s,p) 
where 

(4) 

/    P  y*     P  y  P    y  * 

T    (s,p)    «   1    +  1     [ ±   +  ± +    ..    +  ....         ,f   ,      ,    ]    y      V  (A.  .10) 

2  J  4|  2m?  j  x.=*/2 


r  n3v    (3)  n(2m+l)         (2m+l)  1 

t  (s,p)    *p  +  J  [— 5T—  +  ..  +  —  ]  y  '> 

JI  (2m+l>  2Jx^T/2 


(A. 11) 


The  derivatives  of  y,  and  y  can  be  obtained  with  the  aid  of 

1       2 


2 

■ — ~  +  (b  -  s  cos  x)y  =  0  (A, 12) 

dx^ 

(2m)        (2m+l) 

1  2 

Then  the  ratios  — — » and  --— — —  for  x  -  IT/2   are  know?',  functions  of 

yi  Y2 

and  S  and  since  the  asymptotic  expansion  of  b  is  given  by 


b  m   v  vrs"  .  2lLt_i  -   y3   +  3^   __  _^v_4_+,34V2__+_9    ^  ,A   13 

8  26>/s"  210s 

T  (S,P)  and  T  (s;p)  can  be  determined.   Therefore  y  (s^  ff/2)  can  be  calculated 

IT  7T 

from  the  above  expressions,  only  if  y  (s  ■-  -  P)   y0(s,  —  -  p)  can  be  calcu- 

lated.   Further^  the  computation  of  the  successive  derivatives  become  cumber- 
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some  and  the  accuracy  of  the  results  obtained  is  less  than  for  very  small 
values  of  x. 


(b)   Evaluation  of  N  ; 

n 


By  definition 

2ir 


Nn  =   J   [Sen(s5n)]2  dx 


Evaluation  of  N  presents  the  same  difficulty  as  that  of  Se  (s.x) 

n  nv  >    ■ 

since  the  asymptotic  expression  given  before  by  Equation  (Aw3)  is  not 

valid  over  the  entire  range  of  x  given  by  0  <C  x  ^  2fl\ 

(2) 

(c)   Evaluation  of  He     (s.x): 

(2) 

He   v  •  (s,x>  s  Je    (s,x)  -  j  Ne  (s,x)  (A. 15) 

£  i j  n  n 


The  asymptotic  expansions  of  Je   and  Ne   are  available  and  are  given  by 


Je    (s,x)~[F    (x)    cos   <p  +  F    (x)    sin  <ff\ 


l/4  rr^TT  ^r 

s  vcosn  x 


(A. 16) 


and 


Ne    (s,x)~[F    (x)    sin  (p  -  F    (x)    cos   <£j 
\/cosh  x 


a/4 

s  vcosh 


(A. 17) 


where 


<P  s=    \/s    sinh      x  -    (2n   +   1)    tan    L    (tanh      x)  (A.  18) 
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and 

F()(x)  -  1  +  — r   +  ¥  [ _] 

4  \|  s  cosh  x       512  cosh  x    512  cosh  x 


1    -(v  +14V  +33^)    2^  +124V  +1122^    3^+29^  +1627V, 
2048  cosh  2x      2048  cosh6x      2048  cosh6  x 


(A. 19) 


fi(x)  .  irfo  si*^x  +  ™^inJiiL_   [v3  +  3v  +  4v3  *  f V] 

16  vd   cosh    x         128s   cos  h    x  cosh  x 


sin  hx    r5^4-i-34V2+9         (v6-47v4+667V2+2835)         v64505^4+12139V2+10395-. 
s3/2         2048  cos  h  x  24576   cosh4x  24576    cosh     x 

+    ...  (A. 20) 


In  this  case,  also,  the  accuracy  of  computation  depends  on  the  values  of 

s  v   and  x  as  in  the  case  of  Se  (s,x).   Further,  the  asymptotic  expansion 

n 

for  Ne  (s  x)  is  not  valid  when  x  =  0 

n   ' 

(2) 

The  possible  way  that  the  functions  He     (s,0)  can  be  computed  are 

from  the  expansions  of  Je   and  Ne   in  the  product  of  Bessels  Functions 

n       n 

given  by 

Je2n(s,0)  =  (~l)n  v/jj7i  £   (-1)   De2k  Jk(-~)  ^(""j*)         (A, 21) 

k=0 

Deo 
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k 


Je2n+1(5,0)    =    (-!)"    v/TT/2     £      (-1)      ^e2k+l[2Jk+lC^W^>    Jk<^/2)] 


*=2 (A.22) 

De 
1 


oo  

Ne      (s,  0)    =    (-l)n    sfTft     E   (-1)      De        J    (  /s/2) 
Zn  Zk.     k 

2z2 1 _  (A. 23) 

Deo 


oo  

Ne  (s,0)    =    (-I)"     Vfl/2     E      (-D    ue  [Y,     _  (  /&72)  J,(  Vfe/2)+Y    (  \T&j2)  J       (7 

2n+l  .    n  2k+l      k+1  k  k  k+l 


D81 


(A. 24) 


Calculations  for  the  Image  Pattern  of  a  Smooth  Parabola 

Let  aFb  be  the  section  of  a  smooth  parabolic  reflector  (Figure  32) 
with  0  as  focus  and  f  as  focal  length.  A  plane  wave  polarized  in  the  z 
direction  is  incident  at  an  angle  6  with  the  x  axis.  It  is  required  to 
find  the  electric  field  on  the  line  defined  by  x  =  0. 

Consider  any  point  P(x,y)  on  the  parabola.   The  path  difference  along 
incident  wave  propagation  vector  between  P  and  F  is  given  by 

FQ  =  FS  +  SQ  =  FS  +  QR  =  FQ  cos  6  +  PQ  sin  6 

o 
=  f  tan  4J/2  cos  6  +  p  sin  ^    sin  6 

=  f  tan  4y2  (tan  4^/2  cos  6+2  sin  6)^  since  p  =  f  Secl/2  ,'« 
current  induced  at  P  due  to  the  incident  plane  wave  is  given  by 

T,T3,   A   2      _,.._    <•,   ikf(cos  6  tan2  <\>/2   +  2  sin  6  tan  4>/2) 

I  (P)  =  z   =  cos  (4V2  +  0)  e  .   Einc 

(A.25) 
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FIGURE   32 


At  F    ,    the   field  due   to   P   is   given  by 

i    k*  (2)         / 

Ez(F   )    »  _ I(p)     Ho         [kp   ] 

4iu€ 


4"€  rj 


cos 


2 
■rt"2        S>    ikf^cos    6    tan      ^/2    +    2    sin   6    tan  ^/ 


2) 


7T 


;-i(kp  -   JT/4) 

sfkP 


mc 


(A. 26) 


97 


By  using  asymptotic  expansion  for  H    [hP  ]  and  noting  that  P  ~  P 

and  kP  »  1.   Coordinates  of  P  are  given  by 

x   =  f  tan2  ty/2,    y   =  2f  tan  ^/2 

dxl  2         dyi         2 

— —  =  f  tan  4y2  sec  4,/2  and  — —  =  f  sec  4'/2.   Hence  a 
dy  dx  ' 

small  length  dS  around  P  is  given  by 


dx  2   dy 

dS  =    (—4   +  (-— i)  •  dMJ  =  f  sec  ^/2  d^ 

dx       dx 

Field  at  F'due  to  the  parabola  is  given  by 

91 
(F/u  i  JT/4       "  cos(4V2  +  6)   -ik[p-f  tan^/2  cos6+2  sin6)  ]  •  d^. 


-z 


-e    cos  4»/2 


where 


.2    4 


p'  =    f   sec  4V2  +  y'   -  4f y1  tan  ^/2, 


A  Brief  Discussion  on  Aberrations 

An  ideal  optical  system  has  a  one  to  one  correspondence  and  a  linear 
relationship  between  an  object  and  its  image.   The  rays  from  a  point 
source  will  converge  to  a  point  in  the  image  space  called  its  conjugate 
point  or  Gaussian  image.   In  terms  of  geometrical  optics  the  corresponding 
wavefront  must  be  spherical  in  the  image  space.   Any  defect  that  will 
cause  a  deviation  of  the  wavefront  from  the  spherical  surface  is  called 
the  aberration  of  the  system  and  will  result  in  an  imperfect  image. 

Consider  a  perfect  optical  system  as  shown  in  Figure  33  where  the 
Gaussian  image  of  a  general  object  point  P,  is  P.   Let  P  be  the  origin 
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FIGURE  33 

of  the  cartesian  coordinates  (x^  y}    z)  where  the  z-axis  passes  through  the 
center  0"  of  a  portion  of  a  spherical  surface^  representing  the  wavefront 
at  the  aperture  in  the  image  space^  the  center  of  the  sphere  being  P' 
and  radius  R.   Let  0  P/  the  displacement  of  the  Gaussian  image  from  the 
axis  of  the  system  be  <T  .   Consider  a  point  P  on  the  wavefront  which  is 
uniquely  described  by  two  cartesian  coordinates  x"  and  y"   or  r  and  <j>  of 
the  cylindrical  coordinate  system.   Then 


=  r  cos 


* 


and 


Z_  =  r  sin  § 


(A. 28) 


where  a  is  the  radius  of  the  aperture. 

For  a  perfect  system^  the  wavefront  in  the  image  space  is  given  by 


/'  2     //  2     ;/  2     2 
x    +  y    +  z    =  R 


(A. 29) 


99 

For  a  system  with  aberration,  R  will  no  longer  be  a  constant.   Let  the 
aberration  function,  V((7,  x',  y '') ,  be  the  deviation  of  this  wave  front 
from  a  spherical  one.   It  is  obvious  that  the  deviation  at  P//(x*}    y  ) 
is  a  function  of  x  ,  y".   The  deviation  of  a  given  system  also  depends 
upon  the  displacement  of  the  object  from  the  axis,  or  more  elegantly, 
upon  0,  the  displacement  of  the  corresponding  Gaussian  image. 
Then 

x"2  +  y'2   +  z"2  =  [R  +  V«7,  x",  y")]2  (A. 30) 

2 
Assuming  that  (V/R)   «  1  we  get 

x    +  y"  +    z"2    -  R   -  2  RV(CT,  x",  y;/)  =  0  (A.  31) 

Then  the  ratio  of  the  direction  cosines  of  a  ray  passing  through  the 
point  P  is  given  by 

x  -  R^  :    7»  -   r^:  «"  (A.32) 

3x         9y 

Hence  the  equation  of  the  ray  is  given  by 


x-x';     v-y"  z-z'7 

=  (A. 33) 


Joy  -  ~~W 
x  -R —   y  -R — /    z 

3x       9y 


The  intercept  of  the  ray  with  the  image  plane  defined  by  z  =  0  is 


8v 

3v 

x  =   R  — 

and 

y  - 

=   R  — 

dx' 

V 
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Changing  into  polar  coordinates 

R  ,  .     9v  ,  cos  <♦>  9Vi  ,A  od\ 

=  rlsin<t>  737  +  — 7*  ml  (A*34) 


x 

a 


Br     r  "  3^J 


and 


R  .      8V   sin  4,  8V.  ,A  __, 

y  =  -[cos  +  5--— 2^]  (A'35) 

The  locus  of  the  intercepts  of  the  rays  from  a  zone  r  =  constant  on 
the  wavefront  will  be  a  closed  curve  on  the  plane  z  =  0.   This  curve  is 
called  the  aberration  curve. 

The  aberration  function  V  can  be  expanded  in  a  power  series  of  the 

2   2 

rational  invariants  0"  f    r  and  0"r  cos  <|>  f    each  term  representing  a  single 

aberration.   A  typical  term  in  the  expansion  is  given  by 


/   ff2i+m  rn  cogm 
inm  Y 


or  by 


b.    0"  +m  rn   cos  m  <j>  (A. 37) 

inm  ~ 


where  t,    n  and  m  are  integers  >  0,  while  n  >  m  and  n  -  m  is  even.   b„ 

1  —  '  inm 

is  the  aberration  coefficient.   Following  Nijboer,  we  use  the  second 
expansion  with  the  following  classification  of  aberrations  (lower  order) 
as  shown  in  Table  2. 

Substituting  (A. 37)  into  (A. 34)  and  (A. 35)  we  get  an  aberration 
curve  due  to  a  circular  zone  r  =  const  on  the  wavefront  given  by 

.     rZl+m   R  f     n+m  n-1  ^    n-m  n-1   .   .   _  x  ,  ) 

=  binm  °^    a  1  "  ~2~  *  si*1^"1)  4*  +  "J"  r    sin  (m+1)  4>  f 

(A. 38) 


101 


y  =  b„   (J 

J  inm 


2i+m  R  (  n+m  n-1 


n-m  n-1 


cos(m-l)  (J)  +  —--   r   "  cos(m+l)  fy 


(A. 39) 


when  V  =  0,  corresponding  to  the  perfect  case  all  rays  converge  to  P. 


TABLE  2 


tt=l 

n=2 

n=3 

n=4 

n=5 
n=6 


curvature  of  the 
field  (r2) 


1st  spherical 

4 
aberration  (r  ) 


2nd  spherical 

c 

aberration  (r  ) 


distortion  (r  cos  <fy) 


3 

1st  corns  (r  cos  <j>) 


5 

2nd  coma  (r  cos  §) 


2 
1st  astigmatism  (r  cos  2d>) 


4 
2nd  astigmatism  (r  cos  2§) 


etc 


3rd  astigmatism  (r   cos  2^,) 


An  examination  of  (A. 38)  and  (A. 39)  will  reveal  that  the  main  features  of 
the  aberration  curve  e.g.,  the  general  shape,  symmetry  etc  are  determined 
by  the  value  of  m,  while  n  determines  the  details  of  the  curve  and  the 
particular  value  of  i    seems  to  be  immaterial.   For  example,  when  m  =  0 
(A. 38)  and  (A. 39)  yeild 


2i  R  f   n-1        ') 
X=  binoa   I  lnr    sin<M 


(A. 40) 


and 


.     n2i  R  f   n-1 
y  =  b.    0   -   nr 
Jt  no  a 


cos  (p   I 


(A. 41) 


More  generally  the  order  of  aberrations  are  designated  according  to  the 
value  of  n,  i.e.,  the  power  of  r.   Thus  what  Nijboer  calls  1st  coma  is 
known  as  3rd  order  coma  since  the  value  of  n  is  3  and  so  on. 
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The  aberration  curve  in  this  case  is  a  circle  with  P  as  the  center.   This 
is  called  spherical  aberration.   One  complete  revolution  along  r  =  const 
circle  corresponds  to  one  revolution  along  the  aberration  curve. 
If  m  =  1,  Equations  (A. 38)  and  (A. 39)  reduce  to 


x  =  b„  ,  (T 
inl 


21+1  R 


n-1   n-1 


— —  r    sin  2<f> 


(A. 42) 


^24+1  R  jn+1   n-1   n-1 
y  =  b.  _  (r    —  ]—--  r       +  — —  cos 

inl       a    2         2 


2$ 


(A.43) 


Equations  (A. 42)  and  (A.43)  represent  a  circle  of  radius  proportional 
to  r    in  the  image  plane  but  the  center  of  the  circle  is  displaced  by 
an  amount  also  proportional  to  r    on  the  y  axis.   Hence  the  family  of 
circles  corresponding  to  various  zones  at  the  aperture  has  an  envelope 
which  is  a  pair  of  .straight  lines  as  shown  in  Figure  34. 


FIGURE  34 
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Thus  the  geometrical  optics  image  of  a  point  object  at  P  has  a  comet- 

/ 
shaped  appearance  with  its  Gaussian  image  P  as  its  apex.   It  is  important 

to  note  that  the  aberration  curve  in  this  case  is  not  symmetrical  except 

with  respect  to  the  y-axis.   This  is  called  coma  aberration.   Two  complete 

revolutions  of  the  aberration  curve  corresponds  to  one  revolution  of  r  = 

constant  circle. 

By  putting  m  =  2  in  Equations  (A. 38)  and  (A. 39)  we  obtain  aberration 
curves  which  are  symmetrical  with  respect  to  both  x  and  y  axes.   These 
curves  are  characterized  by  the  fact  that  they  reduce  to  either  a  hori- 
zontal or  a  vertical  focal  line  as  the  receiving  plane  is  moved  on  either 
side  of  the  Gaussian  image.   This  is  called  astigmatism. 

Higher  values  of  m  will  produce  other  aberrations  which  were  first 
studied  by  tfijboer  .   Without  going  into  the  details  of  these  aberrations 
it  can  be  seen  that  the  aberration  curves  will  be  symmetrical  with  respect 
to  both  x  and  y  axes  for  even  values  of  m  and  symmetrical  only  with 
respect  to  the  y  axis  for  odd  values  of  m  . 

The  field  at  any  point  Q  in  the  image  space  for  uniformly  illuminated 
aperture  is  given  by 

rl      r2V        -ik(s+V) 
u(Q)  oc      |     ® r.dr  d$  (A. 44) 


where  s  represents  the  distance  of  a  typical  point  P  on  the  wavefront 
to  Q. 

With  V  =  0,  the  isophoto  lines  on  the  image  plane  are  a  set  of 
concentric  circles  with  P  as  the  center.   When  the  expansion  of  V  i.e. 
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2  b   CT2    rn  cos  m  <j>  is  used  for  V  in  (A.  44)   Nijboer  has  shown  that 
inm 

the  resultant  diffraction  patterns  are  symmetrical  with  respect  to  both 
x  and  y  axes  only  for  the  terms  corresponding  to  both  spherical  aberration 
and  astigmatism  and  higher  aberrations  involving  even  values  of  m. 
Whereas  the  patterns  do  not  possess  the  same  symmetry  for  the  terms 
corresponding  to  coma  and  higher  aberration  terms  involving  odd  values 
of  m.   All  the  terms  having  odd  values  of  m  including  coma  produce 
patterns  symmetrical  with  respect  to  the  y  axis  only. 

For  the  two  dimensional  case,  since  there  is  no  variation  in  the 
<J>  direction,  the  aberration  function  is 

V  =  2  C  yn  (A. 45) 

n 

2i 
where  C   includes  the  factor  CT  .   It  can  be  seen  from  the  foregoing 

considerations  that  all  terms  with  even  values  of  n  represent  spherical 

aberration  while  all  terms  with  odd  values  of  n  represent  coma  aberration. 

Astigmation  is  absent. 

Following  the  treatment  of  the  three  dimensional  case  by  Nijboer   it 

can  be  shown  that  the  diffraction  pattern  will  be  symmetrical  with  respect 

to  the  z  axis  for  a  system  with  only  spherical  aberration.   On  the  other 

hand  with  coma  the  pattern  is  asymmetrical  about  z  axis.   Thus  in  the 

case  of  a  two  dimensional  zoned  mirror  the  two  secondary  maxima  of  the 

intensity  distribution  along  y  axis  on  either  side  of  the  principal  will 

be  equal  for  a  system  with  no  coma.   Any  inequality  of  the  two  can  be 

attributed  to  the  coma  effect  only,  while  for  a  three  dimensional  case 

their  inequality  is  due  to  both  coma  and  other  higher  aberrations 
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corresponding  to  all  odd  values  of  m.   Therefore  the  ratio  of  the  magnitude 

of  the  two  secondary  maxima  will  provide  a  measure  of  the  coma  effect 

for  a  two  dimensional  case.   In  case  this  ratio  is  nearly  unity?  the  change 

in  the  magnitude  of  the  secondary  maxima  is  a  measure  of  the  spherical 

aberration. 


ADDENDUM 

At  the  time  when  this  report  was  in  the  printer' s  hands,  a  report  on 
the  "Experimental  Study  of  a  Diffraction  Reflector"  by  J.  H.  Provencher, 
Air  Force  Cambridge  Research  Center  (April  1959)  was  received.   These 
authors  are  pleased  to  add  that  Provencher'  s  experimental  results  are 
generally  in  quite  good  agreement  with  those  obtained  in  this  report, 
although  not  on  an  exact  quantitative  basis,  for  Provencher'  s  investigation 
is  made  on  a  three-dimensional  model  and  its  focal  length  and  aperture  are 

also  different  from  the  example  computed  in  this  report.   For  a  scan  angle 

o 
equal  to  17.5   the  drop  in  gain  as  reported  by  Provencher  is  about  2  db, 

while  for  the  cylindric  mirror  investigated  here  it  is  considerably  better, 

about  0.5  db.   This  might  be  due  to  less  inherent  aberrations  of  a  cylindric 

system  as  explained  in  this  report.   However,  Provencher' s  beam  width  and 

minor  lobe  level  seem  to  deteriorate  slightly  less  than  those  computed  in 

this  report  up  to  17.5  scan  angle.   This  apparent  contradiction  to  the 

last  statement  might  possibly  be  explained  by  the  fact  that  Provencher  s 

F-number  is  larger  than  the  one  considered  here   and  also  his  reflector  is 

illuminated  by  a  flared  horn,  probably  not  uniformly,  as  witnessed  by  the 

measured  half-power  beamwidth  which  is  greater  than  that  of  a  uniformly 

illuminated  dish  of  the  same  aperture.   It  is  also  of  interest  to  point  out 

in  Provencher'  s  measurements  for  normal  incidence  that  (1)  the  gain  of  zoned 

mirror  is  less  than  that  of  a  corresponding  smooth  paraboloid  and  (2)  the 

minor  lobe  level  of  the  former  is  greater  than  that  of  the  latter.   This 

result  agrees  with  these  authors.   The  second  effect  is  already  explained 

in  this  report  (P.  83),  while,  for  the  first  one,  as  the  computation  (Figure 

18)  indicates,  is  due  to  the  diffraction  effect  by  the  edges  of  all  the 


zones.   It  is  therefore  doubtful  that  the  diffraction  effect  in  this  sense 
plays  an  essential  role  in  the  coma-correction  for  a  zoned  mirror  in  spite 
of  the  fact  that  this  effect  tends  to  make  the  two  secondary  maxima 
slightly  closer  to  equality  (See  Figures  18  -  23). 


ERRATA 


Page  9,  line  15,   delete  P  at  end  of  line 

Page  10,  Equation  (3.8)  should  be   E_ .  _ .  ■-  -  E 

Dili      z. 

inc 


k2 
Page  12,  Equation  (3.13)  should  be 


f 

I.(s)g(s,s')ds  + 


iw      J  1 

strip  i 

Equation  (3.15)  should  be  =  -  E 

inc 

Equation  (3.16)  should  have  small  n  above  the   T,       and   £ 

j~l       k^l 
j^i       k*'j 

Page  13.  Equation  (3.18)  should  be       /     i;/(s)  ..... 

strip  i 

Page  14   Equation  (3  19)  should  be    =  E(s') 

Page  18,  first  paragraph,  curl  E  -  -  B 

Page  21,  Equation  (4„19),  the  second  half  of  the  equation  in  the  numerator 

should  read,       Ho(2\t  )Se  Ol)So  01  ) 

TJ  on  no 

Equation    (4,20),    should    read 


+  1    fy,a)Hv    '(kr)dy 

I  z  o 


o 

2  2 
Page  27.  line  7,  should  be   s  »  k  h  f 

Page  29.  line  6,  should  read,    1,   is  the  current 

hB 

Page  37v  the  first  equation  on  the  page  is  Equation  (5,12) 
Page  44,  lint  2,  should  be   Let  A  =   /2kr  cos  a/2 
Page  45,  line  12,  should  be   H^2,[k|x~x'  jjdx7 


Page  46,  line  2,  should  be     H   '[  k(x'-x)  ]dx' 

line  6,  the  end  of  the  equation  should  be  dt  dx 
Page  47,  line  4,  the  equation  should  be 


,.°°   -ik(h+x)t  i30   -ik(h+x)z 

* =__  dt  =  e-ik(htx)    f  « V-  dz 

J    (l  +  t)r  y/t2-!  J  (2+z)r+W 

1  o 


line  8,  the  denominator  should  read   (2+z 


)r+*  /z~ 


r+4 
line  15,  the  denominator  should  read   (2+z)   z  /z~ 

Page  49,  line  7,  insert   equation  number  (5.39) 

Page  58,   Figure  10,  line  2,   Observation  (Angle  of  Incidence  a  =  0) 

Page  70,   Figure  20  should  read,  Image  Pattern  on  the  Focal  Plane 

Page  74,  line  10,  Figure  26  should  be  Figure  31 

Page  78,  zoned  minor  should  read  zoned  mirror 

Page  79,  the  numbers  on  the  left  side  should  be  0,  -1,  -2,  and  zoned 

minor  should  be  zoned  mirror 

Page  80,  zoned  minor  should  read  zoned  mirror 

Page  82,  zoned  minor  should  read  zoned  mirror  and  the  caption  for  Figure  30 

should  be  Minor  Lobe  Level  of  the  Radiation  Pattern  „ „ „ „ 

Page  87 „,  Number  14  and  15,  Kellar  should  be  spelled  Keller 

Page  88,  Number  30  should  read   "Die  Beugung  und  „„.„"  . 

Page  89,  line  6  should  be   (a)  The  Fourier  expansion  of  Se(s,  x)  i, e 

line  7  should  be     Se  (s,x)  ~ 

n  ? 

Page  91,  lines  3,  8,  and  17  asynptotic  should  be  spelled  asymptotic 

Page  91,  line  9  should  be  with  n  -  0   i,e  V  »  i  . . . . 

Page  92,  Equations  (A. 10)  and  (A„ll)  should  read  respectively 

1A2' 

x=U"/2  /       x=3T/2 


Ti 
Page  92,  line  17  should  read   . ,,f    only  if  y  (s,  —  -  p),  .... 

Je  (s,x)~[F  (x)  cos  cp  +   F  (x)  sin  (p] 

Page  93,  Equation  (A.  16)  — 

s     /cosh  x 

Ne  (s,x)~[F  (x)  sin  (^  -  F  (x)  cos  (p] 

Equation  (A. 17)  — — 

s      /cosh  x 


dX1.2    /yl,2 


Page  97,  line  6,  should  read   dS  =  ,  (- — )   +  (- — )   „  dd»  =  .  .  .  . 

'  \  dx       dx 


line  7,  should  be   Field  at  F  due 


line  8?  should  read 


9 

E    (F   )    =    f    e1   ff/4       f1    cos(4J/2   4-   6)    e-ik[p-f(tan  ty/2 

J  cos     ty/2 

"ei 


Page   98,    Equation    (A. 28)      —  =   r   sin  <j>      and  —  =   r   cos  <j> 
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